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Certain zero-bias anomalies (ZBAs) in the voltage, temperature and magnetic field dependence 
of the conductance GiV, T, H) of quenched Cu point contacts have previously been interpreted to 
be due to non-magnetic 2-channel Kondo (2CK) scattering from near-degenerate atomic two-level 
tunneling systems (Ralph and Buhrman, 1992; Ralph et al. 1994), and hence to represent an 
experimental realization of the non-Fermi-liquid physics of the T = fixed point of the 2-channel 
Kondo model. In this, the first in a series of three papers (I,II,III) devoted to 2-channel Kondo 
physics, we present a comprehensive review of the quenched Cu ZBA experiments and their 2CK 
interpretation, including new results on ZBAs in constrictions made from Ti or from metallic glasses. 
We first review the evidence that the ZBAs are due to electron scattering from stuctural defects that 
are not static, but possess internal dynamics. In order to distinguish between several mechanisms 
proposed to explain the experiments, we then analyze the scaling properties of the conductance 
at low temperature and voltage and extract from the data a universal scaling function V{v). The 
theoretical calculation of the corresponding scaling function within the 2CK model is the subject of 
papers II and III. The main conclusion of our work is that the properties of the ZBAs, and most 
notably their scaling behavior, are in good agreement with the 2CK model and clearly different from 
several other proposed mechanisms. 
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I. INTRODUCTION 

The study of systems of strongly correlated electrons 
that display non-Fermi-liquid behavior has attracted 
widespread interest in recent years, fueled in part by their 
possible relevance to heavy-fermion compounds |§|,§J5) 
and high-Tc superconductivity materials |^,0J8|. On the 
theoretical front, one of the consequences was a renewed 
interest in various multi-channel Kondo models, some of 
which were predicted by Nozieres and Blandin M to con- 
tain non-Fermi-liquid physics. Some of the most recent 
advances were made by Affleck and Ludwig (AL) (see ||l^] 
and references therein), who developed an exact confor- 
mal field theory (CFT) solution for the T = fixed point 
of the multichannel Kondo models. On the experimental 
front, an experiment performed by two of us (RB) p],pl|, 
that investigated certain zero-bias anomalies (ZBAsjin 
the conductance of quenched Copper nanoconstrictions, 
has emerged as a potential experimental realization of 
the 2-channel Kondo (2CK) model and the correspond- 
ing non-Fermi-liquid physics [plT 
icisms of the 2CK interpretation [15 
mechanisms for the ZBAs have been offered |l^Jl^, the 
2CK scenario has recently received important additional 
suppport from experimental results on ZBAs in constric- 
tions made from Titanium p9[ and from metallic glasses 

00 ■ 

In a series of three papers (I, II, III) we shall present 
a detailed analysis of these ZBA experiments and their 
2CK interpretation. The present paper (I) is a compre- 
hensive review of the ZBA experiments that attempts 
to integrate all experimental results on the quenched Cu, 
Ti and metallic glass constrictions into a coherent picture 
(while postponing all formal theoretical developments to 
papers II and III) . Paper II contains a calculation of the 
non-equilibrium conductance through a nanoconstriction 
containing 2CK impurities, which is compared with the 
Cu experiments. Paper III, which is the only paper of 
the three that requires knowledge of AL's conformal field 
theory solution of the 2CK model, describes a bosonic re- 
formulation |2^ of their theory that is considerably sim- 
pler than those used previously and is needed to derive 
certain key technical results used in paper II. 

Let us begin by briefly summarizing the quenched Cu 
ZBA experiments and how they inspired the theoretical 
work presented in papers II and III. 

RB used lithographic techniques to manufacture 
quenched Cu constrictions of diameters as small as 3 nm 
(see Fig. ||), and studied the conductance G{V,T,H) 
through the so-called nanoconstriction (or point contact) 
as a function of voltage {V), temperature (T) and mag- 
netic field (if). Their constrictions were so small that 
they were able to detect electron scattering at the level of 
individual impurities or defects in the constriction. Since 
the energy dependence of the scattering rate can be ex- 
tracted from the voltage dependence of the conductance, 
such an experiment probes the actual electron-impurity 



scattering mechanism. 

For very small eV/kg and T (< 5K), RB observed 
non-ohmic ZBAs in the V- and T-dependence of the con- 
ductance signals of unannealed, ballistic nanoconstric- 
tions. The qualitative features of these anomalies (such 
as their behavior in a magnetic field, under annealing 
and upon the addition of static impurities), which are 
reviewed in detail in the present paper, lead to the pro- 
posal |jl|] that the ZBAs are caused by a special type of 
defect in the nanoconstrictions, namely two- level systems 
(TLSs). This proposal has recently received strong sup- 
port from a number of subsequent, related experiments 
(briefly reviewed in section VII ) by Upadhyay et al. on Ti 



constrictions 
constrictions 



I9| and by Keijsers et al. on metallic-glass 

Mil- 

There are at least two theories for how TLSs can cause 
ZBAs in nanoconstrictions. In the first, based on Zawad- 
owski's non-magnetic Kondo model |23,E4], the interac- 
tion between TLSs and conduction electrons is described, 
at sufficiently low energies, by the 2CK model (reviewed 
in Appendix B of paper II), leading to an energy de- 
pendent scattering rate and hence a ZBA. In the sec- 
ond, Kozub and Kulik's theory of TLS-population spec- 
troscopy p^ , p^ , the ZBA is attributed to a F-induced 
non-equilibrium occupation of the up per a nd lower en- 
ergy states of the TLSs (see Appendix A3 ). 

Though the two theories make quite similar predictions 
for the shape of the ZBA, they make different predictions 
for the F/T-scaling behavior of G{V, T). Whereas Kozub 
and Kulik's theory predicts that G{V,T) does not obey 
any y/T-scaling relation at all, the 2CK scenario predicts 
H that in the regime T <^ T^ and eV <C fceT^ (where 
Tk is the Kondo temperature), the conductance G{V,T) 
should obey a scaling relation of the form 



G{V,T) - GjO^T) 



F{eV/ksT) 



(1) 



where F{x) is a sample-dependent scaling function. 
Moreover, AL's CFT solution of the 2CK problem sug- 
gested that by scaling out non-universal constants, it 
should be possible to extract from F{x) a universal (i.e. 
sample-independent) scaling curve r(x), and that the 
conductance exponent a should have the universal non- 
Fermi-liquid value a = i in striking contrast to the usual 
Fermi-liquid value |25| of a = 2. Since no calculation had 
been provided in Ref. |Q] to support the statement that 
a = i , its status up to now has been that of an informed 
guess rather than a definite prediction, a situation that 
is remedied in papers II and III. 

A detailed scaling analysis |2| showed that the data 
of RB indeed do obey the above scaling relation, with 
a — 0.5 ±0.05. It should be emphasized that the verifica- 
tion of scaling was a very significant experimental result: 
firstly, the scaling relation (^, by combining the V- and 
T-dependence of G{V, T) for arbitrary ratios of V/T^ con- 
tains much more information than statements about the 
separate V- or T-dependence would; and secondly, an 



accurate experimental determination of the scaling ex- 
ponent a is possibly only by a scaling analysis of all the 
data (for a detailed review of this central ingredient of 
the data analysis, see section [Vl|) . Accurate knowledge of 
a is very important, since a succinctly characterizes the 
low-energy critical properties of the physics, enabling one 
to eliminate many otherwise plausible candidate theories 
for the ZBA (such as that by Kozub and Kulik). 

The experimental value for a agrees remarkably well 
with the CFT prediction of a = ^; furthermore, the scal- 
ing curve r(a;) is indeed the same for all three samples 
studied in detail by RB, in accord with the CFT ex- 
pectation that it should be universal and hence sample- 
independent. Thus, this result considerably strenghtened 
the case for the 2CK interpretation of the RB experi- 
ment, within which the experimental demonstration that 
a = 2 i^i remarkably, equivalent to the direct observation 
of non- Fermi-liquid physics. 

Nevertheless, this scaling behavior can conceivably 
also be accounted for by some other theory. Indeed, 
Wingreen, Altshuler and Meir |l^](a) have pointed out 
that an exponent of a = ^ also arises within an alterna- 
tive interpretation of the experiment, based not on 2CK 
physics but the physics of disorder. (We believe that this 
interpretation is in confli ct w ith other important experi- 
mental facts, see section A 1). 

It is therefore desirable to develop additional quanti- 
tative criteria for comparing experiment to the various 
theories. One possible criterion is the scaling function 
r(x). A very stringent quantitative test of any theory for 
the RB experiment would therefore be to calculate the 
universal scaling function r(a;), which should be a finger- 
print of the theory, and compare it to experiment. Papers 
II and III are devoted to this task: r(x) is calculated an- 
alytically within the framework of the 2CK model and its 
exact CFT solution by AL, and the results are compared 
to the RB experiment. When combined with recent nu- 
merical results of Hettler et al. Il2] , agreement with the 
experimental scaling curve is obtained, thus lending fur- 
ther quantitative support to the 2CK interpretation for 
the Cu constrictions. 

The main conclusion of our work is that the 2CK in- 
terpretation can qualitatively and quantitatively account 
for all the scaling properties of the conductance measured 
in the ZBAs of Cu point contacts. The Ti and metal- 
lic glass results add further evidence in support of the 
2CK interpretation, as opposed to other proposed mech- 
anisms. However, we shall note that the 2CK model does 
not account for two phenomena observed in the quenched 
Cu samples. Firstly, the magnetic field dependence of 
the low-bias conductance is rather strong (the 2CK ex- 
planation for the held dependence that was offered in 
Ref. 1^ does not seem to survive closer scrutiny, as dis- 
cussed in Appendix B). Secondly, in many (but not all) 
Cu constrictions the conductance undergoes very sud- 
den transitions at certain voltages Vc |ll| , p6[ (see (Cu.9) 
of section IV) if T and H are sufficiently small. These 
voltages can be rather large (Vc typically ranges between 



5 and 20 mV), implying that some new, large energy 
scale is involved. These two phenomena are not generic 
to TLS-induced ZBAs, however, since they are observed 
neither in metallic-glass constrictions nor in Ti constric- 
tions, which in fact conform in all respects to what is 
expected for 2CK physics. We shall suggest that the 
two phenomena involve (as yet poorly understood) "high- 
energy" physics associated with the strongly-interacting 
system of electrons and atomic tunneling centers. Such 
physics is beyond the scope of the existing 2CK model 
and its CFT treatment, which deals only with the "low- 
energy" aspects of the problem. 

Paper I is organized as follows: In section O we de- 
scribe the fabrication and characterization of nanocon- 
strictions, and summarize some elements of ballistic point 
contact spectroscopy in section III. In section |^ we sum- 
marize the main experimental facts associated with the 
ZBA in the Cu samples, which we state in the form of 
nine properties, (Cu.l) to (Cu.9). The 2CK interpre- 
tation is presented in section M, where its assumptions 
are summarized and critically discussed. Section VI con- 
tains a scaling analysis of the G{V,T) data a.t H = 0. 
The related ZBA experiments on Ti and metallic-glass 
nanoconstrictions are discussed in section VII. Finally, 



we sum mariz e the results and conclusions of this paper in 
section VIII Appendix A describes experimental argu- 



ments for ruling out a number of conceivable explanations 
for the ZBA that could come to mind as possible alter- 
natives to the 2CK scenario. In Appendix B we discuss 
possible sources of magnetic field dependence in the 2CK 
scenario, concluding it is essentially iJ- independent. 



II. THE NANOCONSTRICTION 

A schematic cross-sectional view of a typical nanocon- 
striction (often also called a point contact) is shown in 
Fig. n^. The device is made in a sandwich structure. The 
middle layer is an insulating Si3N4 membrane. This con- 
tains in one spot a bowl-shaped hole, which just breaks 
through the lower edge of the membrane to form a very 
narrow opening, as small as 3 nm in diameter. This open- 
ing is so small that the resistance signal, measured be- 
tween the top and bottom of the structure, is completely 
dominated by the region within a few constriction radii 
of the opening. Hence the resistance is sensitive to scat- 
tering from single defects in the constriction region. 

To obtain the bowl-shaped hole in a Si3N4 membrane, 
electron beam lithography and reactive ion etching are 
used in a technique developed by Ralls [g^l (the details 
relevant to the present experiments are described in Ref. 
jnl, section 2.2). In ultra- high vacuum (< 2x10"^" torr) 
and at room temperature the membrane is then rotated 
to expose both sides while evaporating metal to fill the 
hole (thus forming a metallic channel through the con- 
striction) and coat both sides of the membrane. A layer of 
at least 2000 A of metal (Cu or Ti in the work described 



below) is deposited on both sides of the membrane to 
form clean, continuous films, and then the devices are 
quenched (see property (Cu.l) in section IV). 



III. 



BALLISTIC POINT CONTACT 
SPECTROSCOPY 



A constriction is called ballistic if electrons travel bal- 
listically through it, along semi-classical, straight-line 
paths between collisions with defects or the walls of the 
constriction. This occurs if two conditions are fulfilled: 
Firstly, it must be possible to neglect effects due to the 
diffraction of electron waves, i.e. one needs l/kp <C a, 
where a = constriction radius. Secondly, the constriction 
must be rather clean (as opposed to disordered): an elec- 
tron should just scatter off impurities once or twice while 
traversing the hole. One therefore needs a <^ I, where I 
is the electron mean free path. 

The quenched Cu ZBA-devices of RB reasonably meet 
both conditions: firstly, for Cu l//cp ~ O.lnm, whereas a 
is of order 2-8 nm [as determined from the Sharvin for- 
mula for the conductance, Eq. (^]. Secondly, for clean, 
annealed devices / ^ 200 nm (as determined from the 
residual bulk resistivity). For devices containing struc- 
tural defects, I is reduced to about I ~ 30 nm [see 
(Cu.4)], which is still about twice the constriction di- 
ameter. Thus, we shall henceforth regard the quenched 
Cu ZBA-devices as ballistic constrictions. 

Some aspects of the theory of transport through bal- 
listic constrictions pq , |2^ are reviewed in Appendix A of 
paper II. Here we merely summarize the main conclu- 
sions. 

The differential conductance has the general form 



G(V) 



dI{V) 



dV 



Go + AG{V) 



(2) 



The constant Go, the so-called Sharvin conductance, 
arises from electrons that travel ballistically through the 
hole without scattering. Sharvin showed that for a round 
hole, 



Go = a^e^mep/(27rr). 



(3) 



where a is the radius of the hole, and hence the mea- 
sured value of Go can be used to estimate the size of the 
constriction. 

Any source of scattering in the constriction that 
backscatters electrons and hence prevents them from bal- 
listically traversing the hole gives rise to a backscatter- 
ing correction AG. If the electron scattering rate T^^{e) 
is energy-dependent, AG{V) is voltage dependent. In 
fact, one of the most important characteristics of ballis- 
tic nanoconstrictions is that the energy dependence of 
T~^(e) can be directly extracted from the voltage depen- 
dence of AG{V), which implies that ballistic nanconstric- 
tions can be used to do spectroscopy of electron-defect 
scattering. 



If, for example, the voltage is large enough to excite 
phonons (> 5mV for Cu), the I-V curve is dominated by 
electron-phonon scattering. In this case, it can be shown 
that at T = 0, AG{V) = - {Ae'^m'^VFa^/iTTTi'^) T'^eV), 

where T^^(e') = L d£a^Fp{e) is the relaxation rate 
for an electron at energy e' above the Fermi surface. 
Thus, due to phonon-backscattering processes, the con- 
ductance of any point contact drops markedly at volt- 
ages large enough to excite phonons [T^ > 5 meV for Cu, 
see Fig. p|(a)]. Furthermore, the function a^Fp{eV), the 
so-called point contact phonon spectrum, can be directly 
obtained from dvAG{V). For any clean, ballistic Cu 
nanoconstriction, dvAG{V) should give the same func- 
tion a^Fp(eV), characteristic of the phonon spectrum, 
and indeed nanoconstriction measurements thereof agree 
with other determinations of a^Fp. However, the am- 
plitude of the phonon-induced peaks is reduced dramati- 
cally if there is significant elastic scattering due to defects 
or impurities in the constriction region, as has been mod- 
eled theoretically [Q and demonstrated experimentally 
[ pl| . Therefore, comparing the point contact phonon 
spectrum of a given point contact to the reference spec- 
trum of a clean point contact provides an important and 
reliable tool for determining whether the point contact is 
clean or not. 

For voltages below the phonon threshold {V < 5mV 
for Cu), the T^-dependence of AG{V) is due to scattering 
off defects. For a set of defects at positions Ri, with an 
isotropic, elastic, but energy-dependent scattering rate 
T^^{e), the backscattering conductance has the form |lj] 

/■oo 

AG{V) = -{T{Q)e^/h) dLu[-d^fo{nuj)] (4) 

J — OO 

i 

We factorized out the constant r(0)e^//i to ensure that 
AG has the correct dimensions and order of magnitude. 
We assume that the resistance contribution from each 
defect may be calculated independently - that is, that 
quantum interference for electrons scattering from multi- 
ple defects may be ignored. The Oi and hi are (unknown) 
constants of order unity that characterize all those de- 
tails of scattering by the i-th impurity that are energy- 
independent and of a sample-specific, geometrical nature. 
The hi account for the fact that the probability that an 
electron will or will not traverse the hole after being scat- 
tered off the i-th impurity depends on the position of the 
impurity relative to the hole. The a^ account for the fact 
that impurities that are at different positions Rt in the 
nanoconstriction feel different effective voltages (because 
the amount by which the non-equilibrium electron distri- 
bution function at Ri differs from the equilibrium Fermi 
function fo depends on Ri). 

In spite of the presence of the many unknown con- 
stants Qi, hi, we shall see that it is nevertheless possible 
to extract general properties of t^^{£) from the mea- 



sured AG{V, T) data. For example, from Eq. (^ one can 
deduce that if 

-V rJ.^ -i/n T-^ / In [max(T, e)] , 
T {e,T)-T (O,r)oc| T"ri^ir/T), ^^> 



then AG{V, T) oc | ^""^^ 
where F and F are scahng functions. 



In [max(T, V)] 

F{V/T) , 



IV. EXPERIMENTAL FACTS FOR QUENCHED 
CU SAMPLES 

In this section we summarize the experimental facts 
relevant to the ZBA in quenched Cu samples. Our in- 
terpretation of these facts is postponed to later sections, 
where some of them will be elaborated upon more fully, 
and where most of the figures quoted below can be found. 

The phenomenon to be studied is illustrated by the 
upper differential conductance curve in Fig. H. Its three 
essential features are the following: Firstly, the differen- 
tial conductance shows a drop for \V\ > 5 mV, due to 
the excitation of pho nons, a process which is well under- 
stood (see section III). Secondly, there are sharp voltage- 
symmetric conductance spikes at somewhat larger volt- 
ages (Vc), called conductance transitions in Ref. |l|,|32| , 
because in the DC conductance they show up as down- 
ward steps with increasing V (see figure |l^ below) . Some 
of their complex properties are listed in point (Cu.9) be- 
low. 

Thirdly, the conductance has a voltage-symmetric dip 
near V = 0; this is the so-called zero-bias anomaly 
(ZBA). As a sample is cooled, the temperature at which 
the zero-bias features become measurable varies from 
sample to sample, ranging from 10 K to 100 mK. This 
paper is concerned mainly with the regime V < 5 mV 
dominated by this ZBA. 

The ZBA is a very robust phenomenon. For decades it 
has been observed, though not carefully investigated, in 
mechanical "spear and anvil" point contacts made from 
a variety of materials, see e.g. |3^]. Even the dramatic 
conductance transitions have probably been seen in early 
ZBA experiments p8| ,nthough their presence had not 
been emphasized there.El 

The advent of the mechanically very stable nanocon- 
strictions employed by RB allowed a detailed systematic 
study of the ZBA. Their findings are discussed at length 



^For example, Fig. 3C of p8| shows a d^I/dV^ spectrum 
with sharp signals, more or less symmetric about zero, that 
are consistent with being derivatives of spikes in the dl /dV 
conductance curve. Note that these signals are too sharp to 
be spectroscopic signals smeared by kT, but are indicative of 
abrupt transitions. 



in ||ll[ and ||3^ . We summarize them in the form of 9 im- 
portant properties of the ZBA in quenched Cu nanocon- 
strictions: 

(Cu.l) Quenching: ZBAs and conductance transitions 
[Fig. ^(a)] are found only in quenched Cu samples, 
i.e. samples that are cooled to cryogenic tempera- 
tures within hours after being formed by evapora- 
tion. They are found in about 50% of such samples, 
and in a variety of materials, such as Cu, Al, Ag 
and Pt. Cu was used in the samples discussed be- 
low). 

(Cu.2) Amplitude: Typical values for G{V = 0) vary from 
2000 to 4000 e^/h. The anomaly is only a smaU 
feature on a very big background conductance: its 
amplitude [Gmax — G(V = 0)] varies from sample 
to sample, from a fraction of e'^/h to as large as 
70e'^/h at 100 mK. It's sign is always the same, 
with G{V,T) increasing from G{0,To) as F or T 
are increased. The sample (#1 in Fig. |7|) showing 
best scaling (see (Cu.6) below) had a maximum 
ZBA amplitude of about 20e^/h. 

(Cu.3) Annealing: 

(a) After annealing at room temperature for sev- 
eral days, the ZBA and conductance spikes disap- 
pear, and the conductance curve looks like that of 
a completely clean point contact [see lower curve in 
Fig. 1(a)] . 

(b) Nevertheless, such annealing changes the to- 
tal conductance by not more than 1% or 2% (both 
increases and decreases have been observed), indi- 
cating that the overall structure of the constriction 
does not undergo drastic changes. 

(c) Upon thermal cycling, i.e. brief (several min- 
utes) excursions to room temperature and back, the 
amplitude of the ZBA and the position Vc of the 
conductance transitions change dramatically and 
non-monotonically [see Fig. ||(a)]. The complexity 
of this behavior suggests that the thermal cycling 
is causing changes in the position of defects within 
the constriction, and that the ZBA is very sensitive 
to the precise configuration of the defects. 

(Cu.4) Effect of disorder: 

(a) If static disorder is intentionally introduced into 
a nanoconstriction by adding 1% or more of im- 
purity atoms such as Au to the Cu during evapo- 
ration, the zero-bias conductance dip and conduc- 
tance spikes disappear completely [see Fig. ™b)]. 
Likewise, the signals are absent in samples for 
which water is adsorbed onto the Si3N4 surface be- 
fore metal deposition (the standard sample fabrica- 
tion procedure therefore involves heating the sam- 
ple to ^ 100°C in vacuum, or exposing it for several 
hours to ultraviolet light in vacuum, before the fi- 
nal metal evaporation is done). 

(b) When a strongly disordered region is created 



near the constriction (by electromigration: a high 
bias (100-500 mV) is appUed at low temperatures 
so that Cu atoms are moved around, a method 
controUably demonstrated in |27y34| , |35|] ), the con- 
ductance shows no ZBA either, but instead smaU- 
amphtude, voltage-dependent (but aperiodic) con- 
ductance fluctuations at low voltage [see Fig. ^(c), 
(d)]. That these are characteristic of strongly dis- 
ordered constrictions and can be interpreted as uni- 
versal conductance fluctuations due to quantum in- 
terference, was established in a separate investiga- 
tion Hi, [|l|, chapter 4], |^. 

(Cu.5) Phonon spectrum: For quenched samples, in the 
point contact phonon spectrum the longitudinal 
phonon peak near 28 niV is not well-defined, and 
the total amplitude of the spectrum is smaller by 
about 15% than after annealing. After anneal- 
ing, the longitudinal phonon peak reappears and 
the spectrum corresponds to that of clean ballis- 
tic point contacts. Both these differences indicate 
(see p. 4) that the elastic mean free path I in the 
annealed samples is somewhat longer than in the 
quenched samples. From the phonon spectrum of 
the latter, I can be estimated (see sectio nJll| ) to be 
I ~ 30 nm [for the sample shown in Fig. 0(a)], still 
more than about twice the constriction diameter 
for that device. Note also that the point contact 
phonon spectrum for a quenched device [Fig. ||(b)] 
is qualitatively very different from that of a strongly 
disordered constriction [Fig. @(d)]. These facts, 
viewed in conjunction with (Cu.3b) and (Cu.4c), 
imply that the Cu constrictions displaying ZBAs 
are still rather clean and ballistic. 

(Cu.6) V/T scaling (to be established in detail in sec- 
tion 0): 

(a) At _ff = 0, the conductance obeys the following 
scaling relation if both V < V^ and T < T^, but 
for arbitrary ratio v — eV/kgT: 



G{V,T) - G{0,T) 

'J'a 



F{v) 



(6) 



Relation (^) allows a large number of data curves 
to be collapsed onto a single, sample-dependent 
scaling curve [e.g. see Figs. g(a) and g(b) below]. 
The departure of individual curves from the low-T 
scaling curve in Figs. ||(a) and ||(b) indicates that 
T^ or r has surpassed the crossover scales Vk or 
Tk- From the data, these are related roughly by 
cVk = iksTji, with Tk in the range 3 to 5 K. 

(b) F{v) is a sample-dependent scaling function 
with the properties F{Q) ^ and F{v) oc v" as 
I! ^ (X), and the scaling exponent is found to have 
the value a = 0.5 ± 0.05. 

(c) By scaling out sample-dependent constants, it 



is possible to extract from F{v) a "universal" scal- 
ing function r[v) [shown in Fig. |ll|(b) below]. T{v) 
is universal in the sense that it is indistinguishable 
for all three devices for which a scaling analysis 
was carried out (they are called sample 1,2 and 3 
below) . 

(Cu.7) Logarithms: For T^ or T beyond the cross-over 
scales Vif or T^, G{V,T) deviates markedly from 
the scaling behavior of (Cu.6) and behaves roughly 
logarithmically: For H = Q and fixed, small T, the 
conductance goes like \nV for V > Vk [Fig. ||(a)]; 
similarly, for H = V = and T > Tk, the conduc- 
tance goes like InT [Fig. |(b)]. 

(Cu.8) Magnetic field: 

(a) When a magnetic field (of up to 6 T) is applied, 
the amplitude of the ZBA in Cu devices decreases 
[see Fig. ||(b)]. The change in amplitude can be 
as large as 24 e'^/h if H changes from to 6 T. 
For sufficiently small H (< IT), at fixed T and 
V — 0, the magnetoconductance roughly follows 
G{H,T) ex 1^1 (see Fig. |l| below). However, the 
available data is insufficient to establish linear be- 
havior beyond doubt, and, for example, would also 
be compatible with a |_ff |^/^-dependence. 

(b) The ZBA dip undergoes no Zeeman-splitting 
in _ff, in constrast to the Zeeman splitting that is 
found for devices intentionally doped with magnetic 
impurities such as Mn [see Fig. §(a)]. 

(Cu.9) Conductance transitions: 

(a) Voltage-symmetric conductance transitions 
(spikes in the differential conductance at certain 
"transition voltages" Vc, see Fig. g) occur only in 
quenched point contacts that show ZBAs, but oc- 
cur in at least 80% of these. The spikes disappear 
under annealing, just as the ZBA does (Cu.3a). 

(b) (i) A single sample can show several such con- 
ductance transitions (up to 6 different VcS have 
been observed in a single sample), (ii) If T and 
H are small (say T < IK, H < 0.5T), K is typically 
rather large, with typical values ranging between 5 
and 20 mV, well above the typical voltages associ- 
ated with the ZBA (i.e. Vc > Vk)- The spikes have 
a very complex behavior as a function of temper- 
ature (T) and magnetic field (iJ), including (iii) a 
hysteretic V^-dependence, (iv) a bifurcation of sin- 
gle spikes into two separate ones {Vci,Vc2) when 
B ^0 (Fig. |l|), (v) the iJ-dependent motion of the 
spike positions VcjH) -^ when H becomes suffi- 
ciently large (Figs. R2, 03), and (vi) a very rapid 
narrowing of the peaks with decreasing T. They 
are described at length, from a phenomenological 
point of view, in Ref. |32 



Any theory that purports to explain the ZBA in Cu 
constrictions must be consistent with all of the above ex- 
perimental facts. An extcntion of this list to include the 



results of the recent related ZBA experiments by Upad- 
hyay et al. on Ti constrictions and by Keijsers et al. |21 1 
on metallic-glass constrictions is presented in section VI] . 
In the next section, we shall argue that the 2CK sce- 
nario provides the most plausible interpretation of the 
above experimental facts. A number of alternative ex- 
planations for the ZBA that could come to mind are dis- 
cussed in Appendix A, but all are found to be in conflict 
with some of the above facts. 



and odd linear combinations of the lowest-lying eigen- 
states of each separate well), whose eigenenergies differ 
by A = (A^ -I- A^)^/^. Ultra-fast two-state systems have 
such a large A^; that A too becomes very large, so that 
at low temperatures only the lowest level governs the 
physics. 



1. Slow Fluctuators 



V. THE 2-CHANNEL KONDO (2CK) 
INTERPRETATION 

In this section, we develop the 2CK interpretation of 
the ZBAs in quenched Cu constrictions. It attributes the 
ZBA to the presence in the constriction region of dy- 
namical structural defects, namely TLSs, that interact 
with conduction electrons according to the non-magnetic 
Kondo model, which renormalizes at low energies to the 
non- Fermi-liquid regime of the 2CK model. We begin by 
briefly recalling in section |V|.A some properties of two- 
level systems (or slightly more generally, dynamical two- 
state systems) in metals. Successes and open questions 
of the 2CK scenario are discussed in subsections M.B and 
M.C, respectively, and its key assumptions are listed, in 
the form of a summary, in subsection |V|.D. 



A. Two-State Systems 

A dynamical two-state system (TSS), is an atom or 
group of atoms that can move between two different po- 
sitions inside a material [B8[. In the absence of interac- 
tions, its behavior is governed by a double- well potential, 
generically depicted in Fig. |6|, with asymmetry energy 
A^, tunneling matrix element A^;. The corresponding 
Hamiltonian is 



H' 



TSS 



HA. 



(7) 



where r^ and t^ are Pauli matrices acting in the two- 
by-two Hilbert space spanned by the states \L) and \R), 
describing the fluctuator in the left or right well. 

Depending on the parameters of the potential, the 
atom's motion between the potential wells is classified as 
either slow, fast or ultrafast, with hopping rates t^^ < 
10**s-\ 10**s~i < r-i < IQi^s-i or t'^ > lO^^s-i, re- 
spectively ||39| . Slow two-state systems, called two-state 
fluctuators, have large barriers and neglibibly small A^^, 
and the motion between wells occurs due to thermally ac- 
tivated hopping or incoherent quantum tunneling. Fast 
two-state systems have sufficiently small barriers and suf- 
ficiently large Aa; that coherent tunneling takes place 
back and forth between the wells. Such a system is known 
as a two-level tunneling system (TLS), because its physics 
is usually dominated by its lowest two eigentstates (even 



The fact that two-state systems in metal nanoconstric- 
tions can influence the conductance was demonstrated 
by Ralls and Buhrman |^,^,^ , who observed so-called 
telegraph signals in weH-anncaled devices (at rather high 
temperatures of 20-150K). These are slow, time-resolved 
fluctuations (fluctuation rates of about lO'^s"^) of the 
conductance between two (or sometimes several) discrete 
values, differing by fractions of e^ /h, which can be at- 
tributed to the fluctuations of a slow two-state fluctuator 
in the constriction region. 

Such telegraph signals were also observed by Zimmer- 
man et al. p^,p| , who studied the conductance of poly- 
chrystalline Bi films, a highly disordered material with 
presumably large numbers of two-state systems. They 
were able to measure the parameters of individual slow 
fluctuators directly, finding values for the asymmetry en- 
ergy A^ ranging from as little as 0.08 K to about IK. 
They also demonstrated that in a disordered environ- 
ment the asymmetry energy of a TLS is a random, non- 
monotonic function of the magnetic field, A^ = A^(i/) 
(as predicted earlier in Ref. M]), and hence can be 
"tuned" at will by changing H . The reason is, roughly, 
that Az depends on the difference 6p = pL — pu between 
the local electron densities at the two minima of the TLS 
potential. Due to quantum interference effects that are 
amplified by the presence of disorder, changes in H can 
induce random changes in Sp and hence also in A^. 

Unfortunatly, experiments on slow fluctuators do not 
yield any direct information on the parameters to be ex- 
pected for fast ones, since their parameters fall in differ- 
ent ranges. 



2. Two-Level Systems 

Fast fluctuators or TLSs presumably have the same 
microscopic nature and origin as slow fluctuators, be- 
ing composed of atoms or small groups of atoms which 
move between two metastable configurations, but with 
much lower barriers. Therefore, they anneal away quicker 
than slow fluctuators, which is why they were not seen 
in the above-mentioned Ralls-Buhrman experiments on 
well-annealed samples mU p. 265]. Also, whereas slow 
fluctuators "freeze out" as T is lowered (which is why 
they don't play a role in the ZBA regime of T < 5K), at 
low T fast fluctuators continue to undergo transitions by 
tunneling quantum- mechanically between the wells. 



A fast fluctuator or TLS interacting with conduction 
electrons is usually described by the non-magnetic or or- 
bital Kondo model, studied in great detail by Zawadowski 
and coworkers |^,^ (it is defined and reviewed in more 
detail in Appendices B and C of paper II; for other re- 
views, see IIHQ): 



H^HTss + Y.^ki. 



cp. 



(8) 
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Here cl creates an electron with momentum k and Pauli 

ka 

spin CT. The terms V° and y^r^ describe diagonal scat- 
tering events in which the TLS-atonis do not tunnel be- 
tween wells. The term V^t'-^ describes so-called electron- 
assisted tunneling processes. During these, electron scat- 
tering does lead to tunneling, and hence the associated 
bare matrix elements are much smaller than for diagonal 
scattering: V" /V^ ~ IQ-^. 

Zawadowski and coworkers showed that the electron- 
assisted term V^t^ renormalizes to substantially larger 
values as the temperature is lowered (as does a similar 
V^ry term that is generated under renormalization) . At 
sufficiently low temperatures (where V^ ~ V^ ~ 1^^), the 
non-magnetic Kondo model was shown B5| to be equiva- 
lent to the standard 2-channcl Kondo (2CK) model, with 
an effective interaction of the form 



i^ri/ = «K /de/de'E E 4. (Re 



a, a crcr' 



(9) 

The two positions of the fast fluctuator in the L- and R 
wells correspond to the spin up and down of a magnetic 
impurity (and L-R transitions to impurity spin flips). 
The electrons are labelled by an energy index e, a so- 
called pseudo-spin index a = 1, 2 (corresponding to those 
two combination of angular momentum states about the 
impurity that couple most strongly to the TLS), and the 
Pauli spin index ct =t) i- Evidently, a plays the role of 
the electron's magnetic spin index in the magnetic 2CK 
model, and since the effective interaction is diagonal in 
CT (which has two values), ct is the channel index. 

This (non-magnetic) 2CK model, with strong analogies 
to the magnetic one, yields an electron scattering rate 
T-^{e,T) with the properties p|,|6| 



T-i(e,T)-T-i(0,r) ex 

In [max (T, e)] if 

ri/2f(e/r) if 



(10) 






(The condition A^/Tjc < T is explained in section |VIC] .) 
Hence^ for T > T^ or A^ /T^ <T <^Tk, it yields [via 
Eq. (^] a contribution to the conductance of Sa{T) oc 
InT or T^/^, respectively. The latter is typical for the 



complicated non-Fermi-liquid physics characteristic of 
the 2CK model in the T ^ Tj^ regime. In this respect 
the non-magnetic 2CK model differs in an important way 
from the (1-channel) magnetic Kondo model, for which 
the low-T scaling is of the Fermi liquid form (ex T^). 



B. Successes of the 2CK Interpretation 

We now turn to an interpretation of facts (Cu.l) to 
(Cu.9) in terms of the 2CK scenario |jl^,pi. Our aim here 
is to sketch the physical picture underlying the scenario. 
Those aspects that require detailed analysis, such as the 
scaling behavior (Cu.6) and magnetic field dependence 
(Cu.8), will be discussed more fully in subsequent sec- 
tions. 

Qualitative features: The cooling and annealing prop- 
erties (Cu.l) and (Cu.3) suggest that the ZBAs are due 
to structural defects or disorder that can anneal away 
at high temperatures [although the well-resolved phonon 
spectrum implies that only a small amount of such disor- 
der can be present (Cu.5)]. This conclusion is reinforced 
by the remarkably complex and non-monotonic behavior 
of the ZBA under thermal cycling (Cu.3c), which indi- 
cates that the ZBA probes the detailed configuration of 
individual defects, not just the average behavior of the 
entire constriction region. Subsequent experiments with 
Ti constrictions have shown that the structural disorder 
is located in the "bulk" of the bowl-shaped hole, not on 
its surface, and that it is caused by geometry-induced 
stress occuring in the metal in the bowl-sh aped part of 
the constriction [see (Ti.ld), section VII A | . 

By assuming that the ZBA is due to fast TLSs, i.e. a 
specific type of structural defect, the 2CK scenario ac- 
counts for all of the properties mentioned in the previous 
paragraph. Property (Cu.4a), the disappearance of the 
ZBA upon the addition of 1% Au atoms, can then be 
attributed to the TLSs being pinned by the additional 
static impurities. 

Logarithms and Scaling: Next, we assume that the 
TLS-electron interaction is governed by Zawadowski's 
non-magnetic Kondo model, which renormalizes to the 
2CK model at low energy scales. This explains a num- 
ber of further facts. Firstly, the non-magnetic nature of 
the interaction explains the absence of a Zeeman split- 
ting in a magnetic field (Cu.Sb). Furthermore, the fact 
that the 2CK scattering rate T~^{e,T) has a logarithmic 
form for e > T^(> T) or T > T^(> e) [see Eq. (|l^)] 
accounts, via Eq. (^), for the asymptotic logarithmic V- 
and T dependence (Cu.7) of G{V,T) for V > Vj,{> T) 
or T > Tk{> V). Thus, we identify the experimental 
crossover temperature T^ (— 3 to 5K) of (Cu.6a) with 
the Kondo temperature of the 2CK model. 

Similarly, the scaling form of T~^{e,T) for e,T ■^Tj^ 
[see Eq. (p^] accounts, via Eq. (||), for the observed scal- 
ing behavior (Cu.6) of G{V,T) for V < V^ and T kT^. 
To be more particular, the very occurence of scaling be- 



havior (Cu.6a), and the fact that the experimental scal- 
ing curve T (v) of (Cu.6c) is universal, can be explained 
(see section VI A ) by assuming that the system is in the 
neighborhood of some fixed point. Assuming this to be 
the 2CK non-Fermi-liquid fixed point, the experimentally 
observed scaling regime can be associated with the theo- 
retical expected scaling regime of A^/T^ < T < Tk and 
V < Tk- Moreover, the non-Fernii-liquid value of a = ^ 
that is th en expected for the scaling exponent (see sec- 
tion VIA) agrees precisely with the value observed for a. 



Thus, within the 2CK interpretation, the experimental 
demonstration of a = ^ is equivalent to the direct obser- 
vation of non-Fermi-liquid behavior. Finally, it will be 
shown in paper II that the shape of the universal scal- 
ing curve T[v) is also in quantitative agreement with the 
2CK model. 

Number of TLSs: Each 2CK impurity in the constric^ 
tion can change the conductance by at most 2e^/hE 
Therefore, the sample with the largest ZBA of 70e^/h 
(sample #2 in Fig. ^ would require up to about 40 such 
TLSs in the constriction. However, this is still only a 
relatively small amount of disorder (correspondigig to a 
density of about IQ-^ TLSs per atom 0, p.277]i). The 
sample that showed the best scaling (sample #1 in Fig. ||) 
had a significantly smaller amplitude of ~ 20e^/h, im- 
plying only about 10 active TLSs (that samples with a 
smaller amplitude should show better scaling is to be 
expecte d due to a smaller spread in A's, see (Ti.3) in 
section [VIIA| ). 



C. Open Questions in the 2CK Scenario 

Having discussed the successes of the 2CK scenario, we 
now turn to questions for which the 2CK scenario is un- 
able to offer a detailed explanation, namely the conduc- 
tance transitions (Cu.9), the strong magnetic field depen- 
dence (Cu.8a), and the microscopic nature of the TLSs. 
We shall point out below that (Cu.9) and (Cu.8a) are not 
generic to TLS-induced ZBAs, and speculate that they 



To see this, we note that in the unitarity limit the scat- 



lannel Kondo defect is pro- 
Eq.(2.20)]), and the pl^e 



tering rate of electrons off a k-ch 

portional to ksm^3 (see e.g. [ [47 , 

shift at the intermediate-coupling fixed point is 5 = Tv/2k 

Thus, in the unitarity limit, the contribution to the resistance 

of a fc = 2 Kondo impurity is the same as for fe = 1, namely 

2e^ /h (the 2 comes from Pauli spin). 

'^For example, the QAQ, constriction studied in M has a diam- 
eter of ~ 13 nm [estimated via the Sharvin formula Eq. (m)], 
and there are 10^ Cu atoms within a sphere of this diameter 
about the constriction. Assuming on the order of ~40 active 
TLSs, their density is therefore roughly of order 10~ /atom. 
Although the constriction is believed to be crystalline, not 
glassy, it is worth noting that this density of TLSs is about 
the same as estimates for the total density of TLSs in glassy 
systems. 



are related and must involve some new "high-energy" 
physics, since (Cu.9) occurs at a large voltage Vc- There- 
fore, our lack of understanding of the latter need not 
affect the 2CK interpretation of the low-energy scaling 
behavior (Cu.6). We conclude with some speculations 
about the microscopic nature of the TLSs, and the like- 
lihood that realistic TLSs will have all the properties re- 
quired by the 2CK scenario. 



1. Conductance Transitions 

The fact that conductance transitions occur only in 
samples that have a ZBA (Cu.9a) suggests [Q that these 
are related to the ZBA: if the latter is phcnomenologically 
viewed as the manifestation of some strongly correlated 
state of the system, then conductance transitions corre- 
spond to the sharp, sudden, "switching off" of the cor- 
relations as V becomes too large. For example, in the 
2CK interpretation, interactions of electrons with TLSs 
in the constriction give rise to a strongly correlated non- 
Fermi-liquid state at small T and V. One might speculate 
that if for some reason a large voltage could "freeze" the 
TLSs, i.e. prevent them from tunneling, this would dis- 
rupt the correlations and give rise to a sudden change in 
the DC conductance and hence a spike in the differential 
conductance. 

At present we are not aware of any detailed micro- 
scopic explanation for the conductance transitions. Note, 
though, that they do not occur in all Cu samples show- 
ing ZBAs. Moreover, recent experiments by Upadhyay et 
al. |19(| on Titanium constrictions and by Keijsers et al. 
| pO| on constrictions made from metallic glasses showed 
TLS-induced ZBAs with properties very similar to RB's 
quenched Cu constrictions, but no conductance transi- 
tions at all (see (Ti.6) and (MG.4) in section VIL where 
these experiments are reviewed) . This suggests that con- 
ductance transitions are not a generic ingredient of the 
phenomenology of ZBAs induced by TLSs. Moreover, in 
the quenched Cu samples, provided that H and T are suf- 
ficiently small, the transition voltage Vc at which the first 
conductance transition occurs usually lies well above Tj^ , 
the scale characterizing the extent of the low-energy scal- 
ing regime of the ZBA [see Fig. 0(a)]. (In other words, 
since they don't occur near zero bias, the conductance 
transitions need not be viewed as part of the zero-bias 
anomaly phenomenon at all, if one restricts this term to 
refer only to the low-energy regime.) 

Thus, there seems to be a clear separation of energy 
scales governing the ZBA and the conductance transi- 
tions. The latter must therefore be governed by some 
new large energy scale due to a mechanism not yet under- 
stood. However, due to the separation of energy scales, 
the conductance transitions need not affect our descrip- 
tion of the low-energy scaling regime of the ZBA below 
Tji (which is ^ eVc/fca) in terms of the 2CK model. 



2. Strong Magnetic Field Dependence 

Since the electron-TLS interaction is non-magnetic, 
i.e. not directly affected by a magnetic field, the 2CK 
scenario predicts no, or at best a very weak magnetic 
field dependence for the ZBA. This agrees with the ab- 
sence of a Zeeman splitting of the ZBA for the Cu sam- 
ples (Cu.8b) (which was in fact one of the main reasons 
for the proposal of the non-magnetic 2CK interpretation 
m). However, it leaves the strong magnetic field depen- 
dence (Cu.8a) as a puzzle. (Two indirect mechanism for 
H to couple to a 2CK system, namely via 7J-tuning of 
the asymmetry energy /S.z{H) and via channel symmetry 
breaking, are investigated in Appendix B; they are found 
to be too weak to account for (Cu.8a), contrary to the 
interpretation we had previously offered H.) 

It is therefore very significant that the experiments 
by Upadhyay et al. on Ti constrictions and by Kei- 
jsers et al. on metallic-glass constrictions show ZBAs 
with only a very weak or even no iJ-dependence [see sec- 
tion ^n[ (Ti.5), (MG.3)], in complete accord with 2CK 
expectations. This suggests that, just as the conductance 
transitions, the strong magnetic field dependence (Cu.Sa) 
of the quenched Cu constrictions is not a generic fea- 
ture of TLS-induced ZBAs. Moreover, Fig. 
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suggests 

that in the Cu samples these two properties might be 
linked, because it shows that the strong i?-dependence 
of G{V = 0, H) is related to the fact that the transition 
voltage Vc decreases to as if is increased (Cu.9b,v). 
(In other words, if the strongly correlated state sets in at 
smaller 14 as F is lowered, the voltage-regime < V < Vc 
in which the anomaly can develop is smaller, so that its 
total amplitude is smaller.) 

Since the main difference between the Cu constrictions 
and the Ti and metallic-glass constrictions seems to be 
that the former contain TLSs with very small A's (see the 
next subsection), whereas in the latter, being disordered 
materials, there will certainly be a broad distribution of 
splittings, we speculate that the conductance spikes and 
strong if-dependence might both be a consequence of 
the very small As occuring in the Cu samples, perhaps 
due to interactions between several TLSs with very small 
splittings. 

Thus, we conclude that attempts (such as those in ||]) 
to explain the -ff-dependence of the ZBA (even at F = 0) 
purely in terms of the 2CK model, which captures only 
the physics at low energies below T^, are misdirected, 
because the 7f-dependence would arise, via the conduc- 
tance transitions, from the "high-energy" physics associ- 
ated with the large scale Vc. 

This interpretation, according to which a magnetic 
field does not directly affect the low-energy physics of 
the phenomenon (only indirectly via its effect on Vc), 
can be checked by doing a V/T scaling analysis at fixed 
but small, non-zero magnetic field. If H is sufficiently 
small that the conductance transitions still occur at rela- 
tively high voltages (i.e. Vc > T^), the scaling properties 



of (Cu.6) should not be affected by having H ^ Q). The 
presently available data is unfortunately insufficient to 
test this prediction. 

The conclusions of this and the previous subsection 
are summ arized in assumptions (A3) and (A4) in subsec- 



tion 



VD 



3. Microscopic Nature of the TLS 

Finally, the 2CK interpretation is of course unable to 
answer the question: What is the microscopic nature of 
the presumed TLSs? Now, ignorance of microscopic de- 
tails does not affect our explanation for why the scaling 
properties (Cu.6) of the ZBA seem to be universal: be- 
cause the latter are presumably governed by the fixed 
point of the 2CK model, any system that is somewhere 
in the vicinity of this fixed point will fiow towards it as 
the temperature is lowered (provided that relevant per- 
turbations are sufficiently small) and hence exhibit the 
same universal behavior, irrespective of its detailed bare 
parameters. 

However, the quality of the scaling behavior implies 
some rather stringent restrictions on the allowed proper- 
ties of the presumed TLSs, because we need to assume 
that all active TLSs (e.g. about 10 for sample #1, which 
shows the best scaling) are close enough in parameter 
space to the non-Fermi-liquid fixed point to show pure 
scaling. 

This implies, firstly, that interactions between TLSs 
(which are known to exist in general p4| , mediated by 
strain fields and changes in electron density), must be 
negligible, because they would drive the system away 
from the 2CK non- Fermi-liquid fixed point. Secondly, 
the fact that scaling is only expected in the regime 
A'^/Tk <T <Tk can be used to estimate that T^ ~ 3 to 
5K and A ~ IK (see section VI for details). Kondo tem- 
peratures in the range of 1-10 K are in good agreement 
with the most recent theoretical estimates for TLSs [Q . 
However, the condition A ~ IK implies that for active 
TLSs the distribution of energy splittings, P(A), must 
be peaked below A :S IK. Since A = (A^ + A^)i/^ both 
the asymmetry energy A^ and tunneling rate A^; must 
be ^IK, a value so small that it needs further comment. 

First note that it is not immediately obvious that val- 
ues of the bare tunneling rate A^; exist at all that allow 
2CK physics: For transitions to be able to take place, 
the barrier between the wells must be sufficiently small, 
but a small barrier is usually associated with a large bare 
Ax , implying a large bare A (and A sets the energy scale 
at which the renormalization flow toward the non- Fermi- 
liquid fixed point is cut off) . Now, for a TLS in a metal, 
the physics of screening can reduce the direct tunnel- 
ing rate A^^ by as much as three orders of magnitude 
under renormalization to T ^ T^ pO| (when tunneling 
between the wells, the tunneling center has to drag along 
its screening clowd, which becomes increasingly difficult. 
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due to the orthogonality catastrophy, at lower temper- 
atures). Thus, the renormalized direct tunneling rate 
can always be assumed very small. Though this implies 
a large effective barrier, it does not necessarily prohibit 
2CK physics: Zarand and Zawadowski WO 51 1 have shown 
2CK physics can be obtained even if A^; = 0, provided 
that the model contains some other channel for inter-well 
transitions, such as electron-assisted transitions via more 
highly excited TLS states (see Appendix C.5 of II). 

More serious is the assumption that the renormalized 
asymmetry energy A^ also be SIK. This may seem very 
small when recalling that in glassy materials, the dis- 
tribution -P(A) for the asymmetry energy is rather flat, 
with A varying over many (often tens of) Kelvins. Note, 
though, that far from being glassy, the constrictions are 
believed to be rather clean (Cu.5), containing almost per- 
fectly crystalline Cu. Therefore, our intuition about the 
properties of TLSs in glasses can not be applied to the 
present system. For example, the TLSs could possibly be 
dislocation kinks. (This would naturally account for the 
disappearance of the ZBA when static disorder is added 
(Cu.4a), since dislocation kinks can be pinned by other 
defects.) Since the dislocation kink would find itself in a 
rather crystalline material, some lattice symmetry could 
guarantee then that the two wells of the TLS are (nearly) 
degenerate and hence assure a small A^ and hence small 
A, etc. 

Moreover, some role might be played by the mecha- 
nism of "autoselection" : This assumes that a given TLS 
will only be "active" , in the sense of contributing to the 
non-trivial V- and T-dependence of the conductance, if 
its (renormalized) parameters happen to be in the ap- 
propriate non-Fermi-liquid regime; if they are not, the 
TLS would only be an "inactive spectator" that only af- 
fects the V- and T-independent background conductance 
Go- Moreover, provided that the distribution ^(A^) is 
not zero near A2 = (which seems very unlikely), there 
should always exist a few TLSs with A ~ IK, since A^, is 
strongly reduced by screening. 

Despite the above arguments, the fact that the effects 
of non-zero As did not show up in the quenched Cu sam- 
ples, necessitating the assumption that all active TLSs 
must have A ~ IK, remains probably the weakest point 
the 2CK interpretation of these samples. Therefore it is 
very significant that the effects of non-zero As were re- 
cently explicitly demonstrated in the Ti constrictions of 
Upadhyay et al: in some samples scaling breaks down be- 
low an energy scale associated with A, which was found 
to be tunable by elect romig r at ion and the application of 
a magnetic field [see section |VII Aj (Ti.4), (Ti.5)]. Thus, 
the new Ti experiments significantly bolster the 2CK sce- 
nario at its hitherto weakest point. 

To shed further light on the effects of non-zero A, it 
would be interesting if ZBA experiments with specific 
type of defects with known parameters could be per- 
formed. 

It should be mentioned that Wingreen, Altshuler and 
Meir (WAM) have recently claimed |5|(a) that the 2CK 



interpretation is internally inconsistent if one takes into 
account the effect of static disorder: using values for 
the coupling constants deduced from the observed Kondo 
temperatures, they concluded that renormalized energy 
splitting A would be dramatically increased (to typi- 
cal values of about lOOK), and in particular that there 
would be zero probability for zero splitting (P(0) — 0). 
However, since their arguments neglected the physics of 
screening (i.e. the strong reduction of A^, under renor- 
malization to lower temperatures), we believe that their 
conclusions, in particular that P(0) = 0, are not persua- 
sive [^ (b) , Q . A critical discussion of their arguments 
is given in Appendix D of paper II. 



D. Summary of Assumptions of 2CK Scenario 

The 2CK interpretation of the ZBA in quenched Cu 
constrictions developed in this section can be summa- 
rized in the following assumptions: 



(Al) The ZBA is due to the presence in the constriction 
region of structural defects, namely TLSs, that in- 
teract with conduction electrons according to the 
non-magnetic Kondo model, which renormalizes at 
low energies to the non-Fermi-liquid regime of the 
2CK model. 

(A2) The TLSs may occur with a distribution of (dressed 
or renormalized) parameters, but all "active" TLSs, 
i.e. those which contribute measurably to the volt- 
age dependence of the conductance in the Cu point 
contacts, must have parameters which cause their 
behaviors to be governed by the physics of the non- 
Fermi-liquid fixed point of the 2CK model. This 
requires that all of these active TLSs do not inter- 
act with each other. Moreover, pure V/T scaling 
will only occur in the window t^ jT^^ < T < T^ 
(in quenched Cu samples it was found for T > 0.4K 
and Tjf > 3 to 5K, implying A :S IK). Note, though, 
that the absence of scaling for T < A^/T^ does not 
necessarily imply the absence of non-Fermi-liquid 
physics, which can still show up as V^'"^ behavior 
for A2/r^ <V <T^. 

(A3) The large magnetic field dependence (Cu.8a) and 
the conductance transitions (Cu.9) of the quenched 
Cu ZBAs are related, but not generic to the ZBA. 
They cannot be explained by 2CK physics alone, 
but involve some new energy scale on the order of 

(A4) A magnetic field does not directly infiuence the low- 
energy physics (y < Vk, T < T^) of the ZBA. 
Therefore, 2CK physics can account for the behav- 
ior of the ZBA in quenched Cu samples at fixed 
H, provided that H is sufficiently small (~ 1 T) 
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that the conductance transitions do not influence 
the 2CK scahng regime (i.e. Vc > Vk)- 



VI. SCALING ANALYSIS OF G{V,T) 

In this section, we present a detailed scaling analysis 
of the data and establish the scaling properties of (Cu.6). 
This is a very important part of the analysis of the ex- 
periment, since the scaling properties were used above 
to eliminate quite a number of otherwise plausible can- 
didate explanations of the ZBA. Nevertheless, since the 
upshot of this section is contained completely in prop- 
erty (Cu.6), this section can by skipped by readers not 
interested in the details of the scaling analysis. 

Of course, the scaling properties (Cu.6) we are about 
to establish are simply experimental facts, independent of 
any theoretical interpretation. Nevertheless, during the 
writing of Ref. |2], these properties were predicted (be- 
fore their experimental verification) on the basis of the 
CFT solution of the 2CK model, and we shall present 
our analysis within this framework. We begin by giv- 
ing in section 



Vl.A the general scaling argument first 



reported in [^ to motivate the scaling Ansatz for the 
conductance G{V,T), and a back-of-the-envelope calcu- 
lation of the scaling function r(a;). A more careful calcu- 
lation, tailored to the nanoconstriction geometry and its 
non-equilibrium peculiarities, is reserved for paper II. 



A. The Scaling Ansatz 

The 2CK model is known M to flow to a non-trivial, 
non- Fermi- liquid fixed point at T = 0, at which the model 
has been solved exactly by Affleck and Ludwig (AL), us- 
ing CFT [|lO[ . This fixed point governs the physics in the 
non-Fermi-liquid regime, namely A^/T^ < T < T^ and 
V <Vk- We shall now show by a general scaling argu- 
ment [pj that the assumption of proximity to this fixed 
point (A2), implies the scaling properties (Cu.6) of the 
conductance G{V,T). 



1. General Scaling Argument 

Consider first the conductance signal Gi{V, T) due to a 
single TLS (labeled by the index i) with T < TJ^ , eT^ < 
kgT^, Ai = 0, but arbitrary ratio eV/kgT. According 
to the general theory of critical phenomena, one expects 
that physical quantities will obey scaling relations in the 
neighborhood of any fixed point. For the conductance in 
the present case, a natural scaling Ansatz is: 



G,(i/,r) = G„-fS,r"r 



A,eV 

{ksTr/p 



fll) 



The parameters Ai and Bi are non-universal, positive 
constants, analogous to the of and hi of Eq. (H), which 



may vary, for instance, as a function of the position of the 
TLSs within the constriction. However, the function r(i;) 
should be a universal function, a fingerprint of the 2CK 
model that is the same for any microscopic realization 
thereof. It must have the asymptotic form V{v) ex v^ 
as ti —^ oo, since G{V,T) must be independent of T for 
eV 3> fcsT. Due to the arbitrariness of Ai and Bi, we 
are free to use the normalization conventions that 



r(o) ^ 1 , v{v) 



^vP 



constant as v 



P 



(12) 



Now, ifV is small enough, its only effect will be to cre- 
ate a slightly non-equilibrium electron distribution in the 
leads. In particular, effects that directly affect the impu- 
rity itself, like ^-dependent strains, or the "polarization" 
of the TLS in one well due to the non-equilibrium electron 
distribution, etc. can then be neglected. In this case, 
which we shall call the weakly non- equilibrium regime, V 
only enters in the Fermi functions of the leads, in the 
form [e/^f^"^^/^) -|- 1]~^, i.e. in the combination eV/ksT, 
implying a = (3. 

For a constriction with several defects, the conductance 
signal will be additivep i.e. (now using a— (3): 



G{V,T)^Go + T'^J2^^'^ 



AeV 
knT 



(13) 



Subtracting G(0,T) from this to eliminate Go then im- 
mediately results in the scaling relation of (Cu.6a): 



The scaling function F{v) defined on the right-hand side 
is non- universal, since it depends on the Ai and Bi. 

This is as far as general scaling arguments will take us; 
a specific theory is needed to predict a. To this end, we 
argue by analogy with the conductivity of a bulk metal 
containing 2CK impurities. There the bulk conductivity 
CT(r) is determined, via the Kubo formula. 



a{T) 



'3to2 J (27r)3 



[-d,Jo{e)]p'T{ej,), (15) 



by the elastic scattering life-time T~^{e) = — 2Im S'^(e) , 
where e = £p—£p, and S^(e) is the retarded electron self- 
energy. I]^(e) has been calculated exactly by Affleck and 
Ludwig, using CFT pa, for the bulk fc-channel Kondo 
problem in the neighborhood of its T = fixed point 
(i.e. for T <C T^). They found that in general r^^ has 
the following scaling form (motivated in paper III, or [[46| , 
P, chapter 8]): 



*To be more precise: the contributions of the impurities to 
the resistance R = G~^ are additive, but since 7? = _Ro -f 
Y,, SR,{V,T), with Ro > SR^{V,T), the form (|l|) follows. 
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T-i(e,{A„J) = -2ImE^(£,{A„,}) 



(16) 



1 + E^'^ 



T"'" tra{e/T) 



The sum on m is over all perturbations to the fixed 
point action that one wants to consider. Each such per- 
turbation is characterized by a non-universal parame- 
ter Am. which measures its strength (and has dimensions 
of E^""% where E„i is the energy scale characterizing 
this perturbation) , a universal scaling dimension a„i and 
a (dimensionless) universal scaling function T^ix). In 
principle all the am and Tmix) (but not the non- universal 
Am) can be calculated exactly from CFT, provided all the 
AmT'"™ are small enough that one is in the close vicin- 
ity of the fixed point. Perturbations with am < or 
> are relevant or irrelevant, respectively, because they 
grow or decrease as the temperature is lowered at fixed 
Am- For all perturbations of interest in this paper, the 
scaling functions have the properties rm(a;) = Tm{—x) 
and Tmix) oc x""* as a: — > oo (the latter property follows 
because the perturbation must become T-independcnt in 
the limit e »T). 

AL have calculated in detail the leading irrelevant cor- 
rection to T^^ for the fc-channel Kondo problem, for the 
case that no relevant perturbations are present. In other 
words they take Am = for all m for which am < 0, 
and consider only the correction corresponding to the 
smallest am > 0, say ai. When referring only to this 
correction, we shall drop the subscript m = 1 and de- 
note the corresponding parameters by Ai = A, ai = a 
and ri(a;) = r(x). They showed that for the fc-channel 
Kondo problem a = 2T^j A = AT~" > (where A is 
a dimensionless number of order unity) and r(a::) < O.eI 
(For an explicit expression for T{x), see M, eq.(3.50)] or 
paper II). 

Considering only this leading irrelevant perturbation, 
it follows immediately from the Kubo formula that for 
the 2CK problem (k — 2, hence a = ^), the bulk con- 
ductivity has the form 



a{T) 



T 

t7< 



1/2 



CTl , 



(17) 



with (Ti > 0. The unusual power law T^ " is a signature 
of the non- Fermi-liquid nature of the T = fixed point. 
For a Fermi liquid, one would have had T^. 

Although the form (p6|) for T~^(e) was derived by AL 
only for a bulk geometry, it is natural to assume that 
it also governs the conductance in the nanoconstriction 



geometry (paper II is devoted to a careful justification 
of this assumption). This implies that the exponent in 
Eq. (^l|) should also be a = i, which completes our 
general-principles motivation for the scaling Ansatz. 



2. Back- of -the- envelope calculation of r{v) 

If one is willing to gloss over important subtleties, it is 
possible to obtain by a simple back-of-the-envelope calcu- 
lation a quantitative expression for the scaling function 
that agrees with that found by more careful means in 
paper II. 

Our starting point is Eq. (Q), which gives the change 
in conductance due to back-scattering off defects in a 
nanoconstriction in terms of the scattering rate t^^{s). 
Now, the main difference between a bulk metal and a 
nanoconstriction is that the latter represents a decidedly 
non-equilibrium situation. However, in the weakly non- 
equilibrium regime, i.e. if the voltage is small enough 
{V < Vk), it is a reasonable guess (which is substantiated 
in II) that the form of the scattering rate of electrons off 
a TLS in the nanoconstriction is not all that different 
when the TLSs are in the bulk. Hence, let us boldly usel 
the equilibrium form for t~^ , namely Eq. (|l6|), in Eq. (jj) 
for AG, thus obtaining (to lowest order in Am): 

G{V,T,{Xm}) ^Go- eyhY,^rnT"'-Jde[-d,fo{e)] 

m 

X ^64[fm(£-iel/a+)-f f(e + ieya-)] (18) 

i 

Now write e/fc^r = x, eV/ksT = v, /„(«) = [e" + 1]^^, 
and define (universal) functions Tmiv) by: 

Jb.mTmha.rnV) = dx[dx fo{x)]tm{x + V /2) . (19) 

Here ^a,m and ^b,m are universal constants, chosen such 
that Tm{v) is normalized as in Eq. (I2h, with rm(oo) > 0: 

rm(0) = 1 , Tm{v) = v"'" + const., as w"" -^ oo . 

(20) 

Using the property rm(a;) = rm(— 2:) in the first term of 
Eq. (Eq), we find 

G{V, T, {Am}) =Go + eV/i^ AmT"'"7,,m 

m 
x'^bi ^ Tm{al'Ja.,rnV)+rm{a^'-fa,7nV) ■ (21) 



^The sign of A is a priori undetermined in the CFT approach 
(see ma, after eq. (3.64)]); however, to conform to the expec- 
tation that the Kondo scattering rate increases as e or T are 
decreased, we need A > 0, since r(a;) < 0. 



®The justification for this assumption is explained in sec- 
tion v.. 2 of paper II; essentially, we assume that the leading 
non-equilibrium corrections to r~^ are of order V/Tk, which 
are negligible in the weakly non-equilibrium or scaling regime. 
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Let us now specialize again to the leading irrelevant 
perturbation, for which only Ai = A 7^ 0. (Since for this 
case T{x) < 0, 7a, 1 = ja and 7b, 1 = 7b can both be chosen 
positive.) In this case, Eq. ( |2l| ) is precisely of the form of 
the scaling Ansatz Eq. (|3|), with a = ^. Thus we have 
found a "derivation" for the scaling Ansatz E Moreover, 
this little calculation has furnished us with an expression, 
namely Eq. (|l9|), for the universal scaling function r{v) 
in terms of the exactly known universal function T{x). 

This, in a nutshell, is all there is to the scaling predic- 
tion. Of course, to back up this result by a respectable 
calculation, considerably more care is required and sev- 
eral technical and conceptual hurdles have to be over- 
come. These are addressed in paper II. And finally, in 
paper III the origin of the exponent a — ^ is explained; 
AL originally derived it using the full machinery of CFT, 
but it in III it is found much more simply and directly 
using a recent reformulation ||22| of their theory in terms 
of free boson fields. 



B. Scaling Analysis of Experimental Data 

In this section we summarize the results of a careful 
scaling analysis of the experimental data yj, based on 
section 6.4.2 of 0, and estabhsh that G(V,T) has the 
properies summarized in (Cu.6) of section IV. 



1. First Test ofT'^^'^ and V^^'^ Behavior 

As a first test of the scaling relation Eq. ( J13| ) , one can 
consider it in the asymptotic limits v —>■ and 00 , in 
which the conductance becomes [using Eq. (p^)]: 

G(0,T) = Go + r"Bs, (22) 

G{V, To) = const + v'^Fo at fixed To < eV/k^ , (23) 

where B^ = J2^B^ and Fo = E^-S^Af. Figs. 0(a) 
and |(b) confirm that G(0^) and G{V, 0) roughly con- 
form to Eqs. (E2h and (23), respectively, with a = i. 
Values for Bs and Fq can be obtained from straight-line 
fits to these data, and are listed in table |. However, the 
quality of these data is not good enough to rule out other 
values of a, ranging from 0.25 to 0.75. 



'^Note though that Eg. (IlII) is actually a little too simplis- 
tic, since in Eq. (Ell) eacn defect gives rise to two terms 
with different af and a~ . Note also that bi, 70 are by 
definition all positive constants, and fYa;) > 0. With the 
choice A > 0, discussed in footnote m, we therefore have 
G(V, T, A) - G(0, T, A) > 0, consistent with experiment. 



2. Scaling Collapse 

A much more stringent determination of a can be ob- 
tained from the scaling properties of the combined V 
and T dependence of the G{V,T), which, according to 
the 2CK interpretation, should follow Eq. (|lj). To check 
whether the data obey this relation, the left-hand side 
of Eq. ( |l4| ) should be plotted vs. v". Provided that the 
correct value of a has been chosen, the low-T curves for a 
given sample should all collapse, with no further adjust- 
ment of free parameters, onto the sample-specific scaling 
curve F{v)ys. v", which should be linear for large v" 
[by Eq. (P3|)]. By adjusting a to obtain the best possi- 
ble collapse, a can be determined from the data rather 
accurately. The 2CK scenario predicts a — ^■ 

The raw data for the differential conductance G{V, T) 
of sample #1 of Fig. 0is shown in Fig. |q(a), for T ranging 
from 100 mK to 5.7 K. Using a = ^, rescaling as in 
Eq. (O) and plotting the left-hand side vs. f^^^, these 
data have the form shown in Fig. ^(b). The data at low 
V and low T collapse remarkably well onto one curve, 
which we shall call the scaling curve. Furthermore, F(v) 
vs. v^'"^ has linear asymptote as f — > 00, in agreement 
with Eq. (H). 

Most of the individual curves deviate from the scaling 
curve for V larger than a typical scale eV^ — ImeV (this 
is also roughly the voltage at which lower-lying curves in 
Fig. ||(a) begin to fall on top of each other) . Likewise, 
the lowest curves in the figure, which correspond to the 
highest temperatures, deviate from the scaling curve for 
almost all V . These deviations from scaling at high V 
and T are expected, since if either V > Vk ot:T > T^, the 
scaling Ansatz is expected to break down. We estimate 
Tji as that T for which the rescaled data already deviate 
from the scaling curve at eV/ksT < 1. This gives T^ > 
5 K for the defects of sampl e 1, w hich are reasonable 
values, as discussed in section VC3, and establishes the 
empirical relation cVk — 2k gT^- 

A somewhat complimentary estimate for the highest 
Kondo temperatures of the TLSs in these samples comes 
from the temperature at which the zero-bias signals first 
become visible as the samples are cooled. For all 3 sam- 
ples featured here, this value is approximately 10 K. 

The quality of the scaling provides an exacting test of 
the exponent in the scaling Ansatz. Using a = 0.4 or 
0.6 instead of 0.5 in Eq. ( |l4| ) produces a clear worsening 
of the collapse of the data (see Fig. 2(b) of [0). As a 
quantitative measure of the quality of scaling, we define 
the parameter D(a), which is the mean square deviation 
from the average scaling curve F{v) = ;^ X)n=i '^"(^) 
(where n labels the different experimental curves, corre- 
sponding to different temperatures T„), integrated over 
small values of \v\: 

D{a)^^J2 dv[F^{v)-F{v)]\ (24) 
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D{a) — would signify perfect scaling. Taking the 5 
lowest T (< 1.4 K) and Vmax = 8 (these are the data 
which a priori would be expected to be most accurately 
within the scaling regime, since they are closest to the 
T = fixed point), one obtains Fig. 0(a). Evidently the 
best scaling of the data requires a — 0.48 ±0.05 [the esti- 
mated uncertainty of ±0.05 comes from the uncertainty 
in the exact minimum in the curve in Fig. ^(a)]. This is 
in remarkably good agreement with the CFT prediction 
of a = 5. 

We have also tested the more general scaling form of 
Eq. (pTl), and have observed scaling for 0.2 < a < 0.8, 
with (/? — 0.5) ~ {a — 0.5)/2, with best scaling for a — 
0.5 ± 0.05. But as argued earlier on page 12, one expects 
a = (3 on general grounds. 

The scaling Ansatz has also been tested on two other 
Cu samples. The rescaled data for sample 2 [Fig. |lO|(a)] 
collapse well onto a single curve at low V and T, for 
a = 0.52 ± 0.05 [Fig. |(b)] and with T^ > 3.5 K. At 
high V and high T the non-universal conductance spikes 
of (Cu.9) are visible. The data for sample #3 do not 
seem to collapse as well [Fig. |l^(b)] (illustrating how im- 
pressively accurate by comparison the scaling is for sam- 
ples #1 and #2). However, we suggest that this sample 
in fact displays two separate sets of scaling curves (see 
arrows), one for T < 0.4 K and one for 0.6 K < T < 5 K, 
with interpolating curves in between. This could be 
due to defects with a distribution of Tj^'s, some hav- 
ing Tk ~ 0.4 K and others having T^ > 5 K. The sec- 
ond (higher-T) set of curves do not collapse onto each 
other as well as the first, presumably because there is 
still some (approximately logarithmic) contribution from 
the T^ ~ 0.4 K defects. 



3. Universality 

If for any sample all the Ai in Eq. (04) were equal, one 
could directly extract the universal scaling curve from 
the data. The curve obtained by plotting 



G{V,T) - G{0,T) 



versus {AeV/kgT) 



1/2 



(25) 



with A determined by the requirement that the asymp- 
totic slope be equal to 1 [compare Eq. dl2)], would be 
identical to the universal curve (r(a;) — 1) vs. x^/^. Such 
plots are shown in Fig. ^(b). The fact that the scaling 
curves for all three samples are indistinguishable indi- 
cates that the distribution of Ai's in each sample is quite 
narrow and is a measure of the universality of the ob- 
served behavior. 

To make possible quantitative comparisons of the data 
with the CFT prediction of Eq. (jfivj ) , we now proceed to 
extract from the data the value of a universal (sample- 
independent) constant [essentially a Taylor coefficient of 
r(i;)], which is independent of the possible distribution 
of Ai 's and i3, 's. 



Consider the regime v ^ 1. As argued earlier, here 
r(w) ~ v^ , and since (3— a— 1/2, with the normalization 
conventions of Eq. (n3) we can write, asymptotically 



(26) 



r{v)-l = v'/^+Ti + Oiv-'/^), 



where Fi is a universal number. It characterizes how long 
the T{v) curve stays "flat" for small v before bending 
upwards towards its asymtotic t;^'^-behavior. It follows 
from Eq. (ph that 



F{v) ^ v^^^Fo + Fi + 0{v-^^^) , 



(27) 



where Fq = T,^BiAy^ and Fi = TiB^. Values for Fq 
and Fi may be determined from the conductance data by 
plotting F versus v^^'^ and fitting the data for large w^/^ 
to a straight line. For samples 1 and 2 we fit between 
(ey/ZcsT)^' ^ = 2 and 3, and for sample 3 (using only the 
curves below 250 mK) between 2 and 2.5. 

Values for Fq, and Fi are listed in Table ||. The un- 
certainties listed are standard deviations of values deter- 
mined at different T within the scaling regime for each 
sample. From these quantities, an experimental deter- 
mination of the universal number Fi — Fi/B^: can be 
obtained; it is listed in Table |. The values of Fi are 
consistent among all 3 samples, in agreement with our 
expectation that Fi should be a universal number. 

Paper II will be devoted to a calculation of the uni- 
versal scaling curve T(v) and the universal number Fi. 
There it is shown that the quantitative predictions of the 
2CK scenario are indeed consistent with the measured 
scaling curve. 



C. Upper Bound on the Energy Splitting A 

The energy splitting A of a TLS is a relevant pertur- 
bation to the degenerate 2CK fixed point. In the lan- 
guage of the magnetic Kondo problem, it acts like a local 
magnetic field, and hence has scaling dimension — ^ (see 
(3.19) of H, or |3|, section 3.4.1 (e)]). Therefore a non- 
zero A implies that the electron scattering rate T~^(e) 
and hence the conductance G{V, T) will contain correc- 
tion terms (to be labelled by ni — 2), given by Eqs. ( |lq ) 

and (pi|), respectively, with 0:2 = —^ and A2 = A2A/FJ' . 
Here A2 is a dimensionless number of order unity, and 
E2 is an energy that sets the scale at which A becomes 
important. Though this scale is not a priori known in 
AL's CFT treatment, we shall take E2 — T^, since no 
other obvious energy scale suggests itself. 

It follows that such A-dependent corrections, which 
would spoil good V/T scaling, are unimportant only if 
X^T"^ < 1, i.e. as long as 



A < (rr^) = 



(28) 



holds for each active TLS. Note that this inequality al- 
lows A to be somewhat larger than the naive estimate 
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that would follow from A < T. As emphasized by Za- 
wadowski |p3| , this somewhat enlarges the window of pa- 
rameter space in which the 2CK scenario is applicable. 

The above analysis enables us to estimate an upper 
bound on the energy splittings of active TLSs occuring 
in the quenched Cu samples. The data for samples 1 
and 2 show pure T^/'^ behavior at ^ = (i.e. absence of 
A-corrections) for T as small as 0.4 K. Taking r^«5 K, 
Eq. Q implies that for any active TLS, A< 1.4 K. This 
upper bound is rather small, and was discussed at some 
length in section V C 3 . 



VII. RELATED EXPERIMENTS 

In recent years, ZBAs have bee n found in a number of 
different nanoconstriction studies 
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It should 

be appreciated that in each the ZBA could in principle 
be caused by a different mechanism. However, two recent 
experiments have found ZBAs that convincingly seem to 
be of the same type and origin as those in the quenched 
Cu constrictions. We review their properties below in the 
form of a continuation of the list of properties compiled 
in section IV , together with their interpretation in terms 
of the 2CK scenario. 



A. Titanium Nanoconstrictions 

Three of us (ULB) Q have recently studied nancon- 
strictions with the same geometry as the quenched Cu 
constrictions of RB, but with the leads made from Ti- 
tanium (Ti). This is a stressed refractory metal, which 
is both more disordered and in a state of higher tensile 
stress than Cu, and thus is a likely candidate to have 
dynamical defects. The following properties were found: 

(Ti.l) General Properties: 

(a) ZBAs occur in more than 90% of the samples. 

— This high rate of occurence is due to the highly 
stressed nature of refractory ruetals (mean free path 
is estimated to be ^ ~ 10 nm)j3 and the consequent 
abundance of TLSs. 

(b) The typical amplitude is about AG ~ 10e^//i. 

— This indicates that just a few (probably less 
than 5) TLSs are involved, since the 2CK model 
implies a maximum AG of 2e^//i per defect (see 
footnote H). 

(c) The ZBA anneals away at room temperatures 



®Note that the mean free path is not mucli larger (if at 
all) than the typical constriction radius (5- 15 nm). Strictly 
speaking, this means that one is approaching the regime in 
which the theory of diffusive, not ballistic, point contact spec- 
troscopy should be used (section [II), but we shall continue 
to use the latter. 



on a time scale of a few days to a few months. — 
This is significantly longer than in Cu samples, be- 
cause Ti has a higher melting temperature, 
(d) Geometry-induced stress: If a dirty insulating 
substrate is used (e.g. with organic contaminants), 
to which the Ti-film does not stick well, the ZBAs 
were absent in almost all the samples. The ZBA oc- 
curs only if there is good adhesion between metal 
film and substrate. Nevertheless, the ZBA is not 
a surface effect (e.g. due to TLSs on the surface 
caused by a mismatch in lattice constants between 
Ti and the substrate). This was demonstrated as 
follows by studying constrictions made from a com- 
bination of Palladium (Pd) and Ti: Samples made 
purely from Pd almost never had ZBAs; even when 
they did, the ZBA annealed away over a period of a 
few hours, in contrast to the much longer annealing 
times for Ti (Ti.lc). This implies that stresses are 
relieved much quicker in Pd than in Ti. Now, for a 
series of samples, first a layer of Ti of thickness 2, 
5, 10 or 25nm was deposited on the both sides of 
the constriction, and thereafter Pd was used to fill 
up the bowl and form the leads. If the ZBA were 
a surface effect, it should have shown up in these 
samples. However, very few of them had ZBAs, and 
even when they did, the ZBA annealed away over 
a period of a few hours, exactly as for the pure-Pd 
samples. However when the Ti layer's thickness 
was increased above 30nm, which in these sam- 
ples is approximately the cross-over thickness for 
filling the bowl sufficiently to form a continuous 
nanobridge through the constriction, ZBAs started 
to appear with annealing times characteristic of 
Ti. This shows that the ZBA's cause is situated 
in the bulk of the bowl and not on its surface. — 
This revealing investigation implies that the ZBA 
results from geometry-induced stress in the metal. 
This stress is provided by the bowl-like shape of the 
hole (but only when good adhesion to the bowl is 
possible) and anneals away very slowly in Ti, but 
rapidly in Pd. This conclusion strongly supports 
the assumption (Al) of section VD that the ZBA 
is caused by stress-induced structural defects. 

(Ti.2) V/T scaling: Some samples show the same scaling 
behavior, [G{V,T) - G(0,T)]/r" = F{eV/ksT), 
as that observed for Cu samples (Cu.9a). Of those 
samples for which a scaling analysis according to 
section ^ll was done, at least three scale very well 
for a ~ 0.5; scaling plots for one of these (called 
sample 4 here) are shown in Fig. |l5|. 
— This is in excellent agreement with the 2CK 
scenario, given assuinption (A2) that the TLS en- 
ergy splittings A are sufficiently small. The scal- 
ing curve has the same universal shape as for the 
Cu constrictions (compare with Fig. ra), with the 
universal number Fi = -0.81 ± 0.10 [see Eq. ^) 
and Table B, but actually scales better, in that 
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the "bending-down" deviations from scaling occur- 
ing in Fig. ||(b) for Cu at large (F/T)^/^ are ab- 
sent here. This makes it somewhat hard to unam- 
bigously determine the Kondo temperatm'e; esti- 
mates give Tft: ~ 10 - 20K. 

(Ti.3) Premature V^''^ saturation: Other samples, such as 
sample 5 of Fig. nBl show for small T a "premature" 
saturation of the y^'^-behavior, implying a very 
marked breakdown of scaling. This can be quanti- 
fied as follows: The conductance (also for sample 4) 
can be fitted, at fixed T, by the phenonienological 
form 



G{V) = G(0) + a[V^ + (2T,)2]i/4 



(29) 



where T^ characterizes the "saturation energy" at 
which the large-voltage y-'^/^-behavior crosses over 
to the flat low-]/ regime. For the scaling sample 4 
of (Ti.2) one finds T^ —T [this is the reason for in- 
cluding the phenomenological factor of 2 before T]: 
in Eq. (|29|)] , meaning sample 4's saturation is due to 
thermal rounding. In sharp contrast, for sample 5 
one finds that for sufficiently low temperatures, Tx 
(=1.43±0.03K) is much larger than T (by a factor 
of almost 20(!) for the lowest T of 76mK). This im- 
plies "premature (non-thermal) saturation" of V^''^ 
behavior as V is lowered. (No such saturation was 
ever observed in Cu samples down to 50 mK). 
— The 2CK scenario attributes [see assumption 
(A2)] premature y^/^ saturation to the presence 
of some TLS with finite energy splittings A ~ r^,. 
At energies below A, 2CK physics "feezes out", 
and pure V/T scaling behavior is destroyed for 
{TTkf/^ < A [Eq. @)], due to the presence of 
a new energy scale. Note, though, that non- Fermi- 
liquid physics is destroyed only at low energies — 
it becomes observable again at sufficiently high en- 
ergies, as demonstrated by the reemergence of V^^'^ 
behavior above T^, in agreement with the last sen- 
tence of assumption (A2). Importantly, this implies 
that even if a constriction contains a large number 
of impurities with a wide distribution of A's, it can 
show yi/^ behavior at sufficiently large V, though 
its scaling will suffer. (Indeed, for the Cu samples 
the quality of scaling was better (Cu.2) for samples 
with smaller ZBA amplitudes, i.e. fewer TLSs.) 

(Ti.4) Electromigration: 

(a) Electromigration, the application of large volt- 
ages (1^=200 mV, J ~ 10^" A/cm ) for short peri- 
ods of time (10 seconds, several times), can cause 
significant changes in the saturation energy T^. For 
sample 5, Fig. 0(a) shows that T^ changed from 
2.3 to 1.4K as a result of such an electromigration. 

(b) Another device (sample 6) had a more com- 
plicated \ow-V behavior [Fig. |l^(b)], character- 
ized by a sum of two terms of the form (29), 



with two distinct saturation energies T^i and Tx2- 
Upon electromigration, they experienced (oppo- 
site!) changes, from 0.9 to 1.5K and 9.7 to 6.8K, 
respectively. 

— Presumably electromigration, which is known 
to controUably cause defect rearrangement [p7| , p5| , 
modifies the parameters of some TLSs, thereby 
changing their As. (Ti.4b) is direct evidence that 
individual defects are responsible for the ZBA; ev- 
idently, two different TLSs dominate the ZBA in 
this particular case. 

(Ti.5) Magnetic field: 

The magnetic field dependence of G{V^ H) is weak 
and random (see Fig. Hq): when H is changed from 
to 5 T, G(V, H) 

(a) changes by less than 2e^//i, with random sign, 
for V <Tx\ and 

(b) is completely iJ-independent for V > T^- 

— The very weak _ff-dependence is consistent with 
2CK expectations (see Appendix B). It shows that 
the strong _ff-dependence observed for Cu ZBAs 
(Cu.8a) is not a ge neric p roperty of the ZBA, as 
discussed in section VC2 and summarized in as- 
sumptions (A3,A4). For Ti constrictions, the en- 
tire iJ-dependence can be attributed to H-tuning 
of A(i?), i.e. to the fact that in disordered ma- 
terials, the TLS splitting is known [^,|l] to be 
a random function A(i?) of H due to disorder- 



enhanced interference effects (see section V A 1 ) . 
iJ-tuning of A(i?) is consistent with the fact that 
the i/-dependent changes in G have roughly the 
same size (Ti.5a) as those due to changes in A in- 
duced by electromigration (Ti.4). It explains the 
random sign of the magnetoconductance (T1.5a), 
and also explains (Ti.5b), because A(_ff) can only 
affect G{V) below T^. 

(Ti.6) No conductance transitions: The conductance 
transitions (Cu.9) that occured in at least 80% of 
the quenched Cu samples have never been observed 
in any of the Ti samples. — This shows that Cu 
conductance transiti ons (C u.9) are not generic, as 
discussed in section V C 1 and summarized in as- 
sumption (A3). 

The above discussion shows that the Ti nanoconstric- 
tions display all the phenomenology expected from 2CK 
impurities: (i) their amplitudes are sufficiently small to 
be attributed to only a very few TLSs; (ii) some of them 



show good V/T scaling with scaling exponent a 



(in) 



others demonstrably show the effects of finite, tunable A; 
and (iv) they lack the puzzling conductance transitions 
and large magnetic field dependence of the quenched Cu 
samples. Thus, they appear to be "custommade" real- 
izations of 2CK physics in nanconstrictions. 
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B. Mechanical Break Junctions made from Metallic 

Glasses 

Keijsers, Shklyarevskii and van Kempen El] studied 
ZBAs in mechanically controlled break junctions made 
from metallic glasses, which are certain to contain many 
TLSs. They observed the following properties: 

(MG.l) Amplitude and shape: 

The ZBA has qualitatively exactly the same shape 
and sign as that of RB's quenched Cu samples, with 
an amplitude of sometimes more than 100 e^/ft.. — 
The large amplitude is to be expected, since metal- 
lic glasses contain a high concentration of TLSs. 

(MG.2) Slow ZBA fluctuations: 

(a) Remarkably, in some samples the amplitude of 
their ZBA fluctuates between two values G{V) and 
G'{V) (or sometimes several) in a telegraph- noise 
fashion on a time scale of seconds, evidently due 
to the presence of one ( or som etimes several) slow 
fluctuators (see section V A 1 ) in the constriction 
region. 

(b) The amplitude of these telegraph-fluctuations, 
AG{V) = IG - G'\, is of order le'^/h or less, and 
depends on V. It decreases from AG(0) to as 
V increases from to between 5 and 10 mV, see 
Fig. |l|. 

(MG.3) Magnetic field: 

The ZBA shows no 7J-dependence. — This is ex- 
actl y as e xpected in the 2CK scenario (see sec- 
tion V_C2), and confirms the conclusion (Ti.5) that 
the strong 7f-dependence of Cu ZBAs (Cu.8a) is 
not a generic property of the ZBA. 

(MG.4) No conductance transitions: The conductance 
transitions (Cu.9) of the quenched Cu samples 
have never been observed metallic-glass constric- 
tions [£^. — This confirms the conclusion (Ti.6) 
that Cu conductance transitions (Cu.9) are not 
generic. 

The ^-dependence of AG{V) implies that the large 
features of the ZBA and the small amplitude fluctuations 
cannot be unrelated phenomena. For example, it is not 
possible to attribute the overall ZBA to a suppression 
in the density of states due to static disorder (analogous 
to the pr opos al of WAM for quenched Cu samples, see 
Appendix A 1 ) , while assuming the additional small con- 
ductance fluctuations to be caused by an independent 
slow fluctuator. The problem with such a scenario would 
be that the amplitude of the fluctuations, though of the 
right magnitude of < e^//i, would be F-mdependent. 

Keijsers et al. state that the large features of their ZBA 
can be explained by invoking either Zawadowski's non- 
magnetic Kondo model (section V A2 ) o r KK 's theory of 
TLS-population spectroscopy (Appendix [A3| ) to describe 
the interaction of electrons with the fast TLSs in their 



system. They propose that the amplitude fluctuations 
can be explained (in either theory) by assuming that the 
TLS-electron interaction strengths \V^ and V'^ in Eq. (||)] 
of some fast TLSs are modulated between two values^ due 
to short-ranged interactions with a nearby slow two-state 
system, when the latter hops between its two positions. 

Zarand, von Delft and Zawadowski recently pointed 
out [ p6| that the maximum switching amplitudes ob- 
served by KSK are so small [AGmax < le^/h for all 
samples in Ref. [gl[) that they seem to stem from the 
parameter-modulations, induced by a slow fluctuator, of 
only one or two fast TLS, where the parameter expected 
to be most strongly modulated is the TLS asymmetry 
energy, which can be assumed to fluctuate between two 
values, Az <^ A'^. This implies that the experiments of 
KSK constitute the first measurements of the conduc- 
tance contributions of individual fast TLSs, and allows 
an unprecedentedly detailed comparison with theory: by 
calculating G{V,Az) (the contribution to the ZBA due 
to scattering off a TLS with asymmetry energy A^) for 
various asymmetry energies, it should be possible to find 
two values A^ and A^ for which |G(V", A^) - G(y, A^)] 
reproduces the measured AG{V). 

Analysing two samples in detail, Zarand, von Delft and 
Zawadowski showed that AG(T^) could not be fit using 
the TLS-population spectroscopy theory of Kozub and 
Kulik (section AS). However, rather good fits (Fig. U9t) 
were achieved using Zawadowski's non-magnetic Kondo 
model (although the analysis does not completely rule 
out that there can be a small KK contribution [pq]). 

This constitutes possibly the most direct evidence yet 
for the applicability of the 2CK model to TLS-induced 
ZBA's in point contacts. Subsequent measurements of 
the response of such break junctions to RF-irradiation 
pTf support this conclusion (although they, too, do not 
completely rule out a small KK contribution). If both 
the V- and T-dependence of AG were known, a V/T 
scaling analysis [g[ would serve as a further test of the 
non-magnetic Kondo scenario. (At present the break 
junctions are not sufhciently stable against mechanical 
deformations under temperature changes to reliably de- 
termine the T-dependence of AG.) 



VIII. CONCLUSIONS 
A. Summary 

This paper is the first in a series of three (1, II, III) 
devoted to 2-channel Kondo physics. We have reviewed 
in detail the experimental facts pertaining to a possible 
realization of the 2CK model, namely the non-magnetic 
ZBA in quenched Cu nanoconstrictions, and also inte- 
grated into our analysis insights obtained from new ex- 
periments on Ti and metallic-glass constrictions. 

We have summarized the various experimental facts 
for the quenched Cu samples in the form of nine prop- 
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erties, (Cu.l) to (Cu.9) (section IV). Properties (Cu.l- 
5), which are of a mainly qualitative nature and very 
robust, place very strong demands on any candidate ex- 
planations of the ZBA: the zero-bias anomalies disap- 
pear under annealing, and hence must be due to struc- 
tural disorder; they disappear when static disorder is 
intentionally added, and hence cannot be due to static 
disorder - instead they must be due to dynamical im- 
purities; they show no Zeeman splitting in a magnetic 
field (Cu.8b), and hence must be of non-magnetic ori- 
gin. These observations lead to the proposal ||l| that the 
zero-bias anomalies are due to nearly degenerate two- 
level systems, interacting with conduction electrons ac- 
cording to the non-magnetic 2-channel Kondo model of 
Zawadowski |49| , which renormalizes at low energies to 
the non-Fermi-liquid regime of the 2CK model. 

We then presented a quantitative analysis of the V/T 
scaling behavior of the conductance G{V, T) = Go + 
T°'F{eV/kgT), which demonstrates unambigously that 
the scaling exponent has the unusual value of a = i, 
in contrast to the usual Fermi- liquid value of a = 2. 
We argued that this too can naturally be understood 
within the phenomenology of the T = fixed point of 
the 2CK model, within which the experimental verifi- 
cation of a = i constitutes the direct observation of a 
non-Fermi-liquid property of the system. Breakdown of 
scaling for larger T and V values is explained too, since 
for these the system is no longer fine-tuned to be close to 
the T = fixed point, thus spoiling the scaling behavior. 
Estimates of T^ in the range 1-5 K, which is reasonable, 
were obtained, as well as an upper bound for the energy 
splitting of all active TLSs of A S kglK. This bound is 



rather small (and has been criticised, see section VC3), 
but is enforced by the quality of the observed scaling. 

The scaling analysis provides sufficiently detailed in- 
formation about the low-energy physics of the system 
that it enabled us to rule out several other candidate 
mechanisms for explaining ZBAs. (An alternative inter- 
pretation of the scaling properties recently proposed by 
Wingreen, Altshuler an d M eir can be discounted on other 
grounds, see Appendix A 1). 

We then reviewed experiments on Ti and metallic glass 
samples in the form of further properties (Ti.l) to (Ti.6) 
and (MG.l) to (MG.4). They provided further strong 
support for the 2CK interpretation: the Ti data demon- 
strated the destruction of scaling in the presence of a non- 
zero, tunable TLS energy splitting A, and the metallic 
glass data allowed the contribution of a single TLS to be 
measured and compared with theory. 

The 2CK interpretation is sufficiently successful in ac- 
counting for the observed phenomenology of the scaling 
properties (Cu.6), that we believe more quantitative cal- 
culations based on this model to be justified. The remain- 
ing two papers in this series, II and III, are devoted to a 
quantitative calculation of the scaling function ^(v), to 
be compared with the experimental curve in Fig. 0(b). 
The final result is shown in section V of paper II, Fig. 6. 
When our results are combined with recent numerical re- 



sults of Hettler et al. ||12|,^ , quantitative agreement with 
the experimental scaling curve is obtained. 

The main conclusion of this investigation is therefore 
that the 2CK interpretation is in qualitative and quanti- 
tative agreement with the scaling properties of the data. 
It can also account for all other observed properties, with 
only two exceptions for the quenched Cu samples, sum- 
marized below. 



B. Open Questions and Outlook 

There are two experimental observations for the 
quenched Cu samples that do appear to lie beyond the 
present understanding of 2CK physics: the conductance 
transitions (Cu.9) and the apparantly related strong 
magnetic field dependence (Cu.8a). A theory of this phe- 
nomenon would be most welcome. However, these ef- 
fects appear to involve either "high energy" effects or ef- 
fects due to interactions between nearly degenerate TLSs 
which are beyond the scope of the present-day single- 
impurity calculations of the 2CK model, which are ap- 
plicable only to the low T- and l^-regime in the neigh- 
borhood of a T = fixed point. However, these two 
experimental effects should not be considered generic to 
the physics of TLSs in nanoconstrictions: they are absent 
in metallic-glass and Ti constrictions, and particularly 
the latter, which display both scaling and the destruc- 
tion thereof by a non-zero, tunable A, seem to be almost 
"ideal" realizations of 2CK physics. 

On the theoretical side, the 2CK model has recently 
been subjected to renewed scrutiny (catalyzed in part 
by its application to the ZBA and the claim that non- 
Fcrmi- liquid behavior has been observed). The main 
point of contention is whether any realistic TLS would 
ever flow towards the non-Fermi-liquid fixed point of the 
2CK model, because of the inevitable presence of rele- 
vant perturbations that drive the system away from this 
point. 

Wingreen, Altshuler and Meir |15[ have argued that 
static disorder could lead to a significant asymmetry en- 
ergy A between the two states of the TLS (a relevant per- 
turbation). We critically discuss their arguments ||lq|(b), 
pO| in Appendix D of paper II, and judge them not to be 
persuasive. More recently, studying a formulation of the 
model that is slightly different from that introduced by 
Zawadowski, Moustakas and Fisher iQ have discovered 
another relevant operator (which was then interpreted 
by Zawadowski et al. 159] to be due to particle-hole sym- 
metry breaking). However, their conclusions have them- 
selves been questioned in Ref. |p9,pQ| i where the prefactor 
of this new relevant operator was estimated to be neg- 
ligably small, and for other technical reasons, some of 
which are mentioned in Appendix D.2 of paper II. 

Certainly, further theoretical work is needed to fully 
understand the stability, or lack thereof, of the T = 
fixed point of the degenerate 2CK model. Both experi- 
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mental and theoretical work would be welcome to better 
understand the nature of the defects giving rise to ZBAs 
in metal point contacts, and the parameters governing 
these defects. Skepticism of the 2CK interpretation of the 
data is not unwarranted, since this is seemingly an exotic 
effect. However, this model has provided a rather com- 
plete account of the experimental observations (Cu.1-7), 
along with accurate predictions of the scaling properties 
of the conductance signals as a combined function of T 
and V. No other existing model, based on more familiar 
physics, has been able to account for all the data. 
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cal depression in the density of states near the Fermi 
surface of the form SN{e - e^.T) oc -T^/'^F{^^), 
where _F is a scaling function. This in turn would re- 
duce the rate at which electrons incident ballistically into 
the disordered region could traverse the sample. The to- 
tal conductance would hence be reduced by an amount 
AG(F,r) (X dv J^^^^^deSNie-e^^T) = 5NiV,T). Due 
to the scaling form of SN, this argument explains the 
scaling property (Cu.6), and in fact the scaling curve 
F{v) of Eq. (o) that it produces is in quantitative agree- 
ment with that of sample 1 (see [^ Fig.l]). According 
to their estimates, this scenario would require a disor- 
dered region of diameter 50 nm (the size of the bowl), 
a mean free path I = 3 nm and a diffusion constant 
_D = 15 cm^/s, i.e. rather strong disorder. 

The WAM scenario is appealing in that it accounts 
for the unusual T^''^ behavior using well-tested physical 
ideas, without having to evoke any exotic new physics 
(such as 2CK non-Fermi liquid physics). However, it is 
at odds with a number of qualitative (and hence very 
robust) properties of the ZBA |1^, section 6.6.1], |jla, (b)]: 

1. According to (Cu.4a), upon the intentional intro- 
duction of static disorder the ZBAs are not en- 
hanced, as one would have expected in WAM's 
static disorder scenario, but disappear completely, 
in contradiction to the latter. 



APPENDIX A: RULING OUT SOME 
ALTERNATIVE INTERPRETATIONS 

In this appendix we discuss a number of conceivable 
explanations for the ZBA that could come to mind as 
possible alternatives to the 2CK scenario. We argue that 
each is inconsistent with some of the experimental facts 
(Cu.l) to (Cu.9) and hence can be ruled out. Most of this 
material is contained in pi|, p3] and fla]. Nevertheless, 
since some of these arguments have been the subject of 
some controversy [n5l, they deserve to be restated and 
summarized here, for the sake of completeness and con- 
venience. 



1. Static Disorder 

Could the ZBA be due to static disorder? For example, 
one could consider attributing the decreased conductance 
near F = to either weak localization due to disorder ||6l| ] 
or disorder-enhanced electron-electron interactions [|62|| . 
In fact, the latter possibility (first mentioned, but deemed 
implausible, in Ref. ||l2[), was recently advocated ||l^] 
by Wingreen, Altshuler and Meir (WAM) (these authors 
also offer a critique [|5[ of a crucial assumption of the 
2CK scenario, which we discussed in section VC3). 



WAM made the interesting observation that if just 
the region of the device near to the point-contact ori- 
fice were highly disordered, this would give rise to a lo- 



2. The quenched Cu constrictions actually are consid- 
erably cleaner than is assumed in WAM's scenario, 
as can be seen from three separate arguments: 

(a) According to (Cu.5), a direct estimate of the 
mean free path, based on the point contact phonon 
spectrum (a reliable and well-tested diagnostic 
method llgj^^) suggests I ~ 30 nm instead of 
WAM's 3 nm. 

(b) WAM attempted to explain (Cu.3a), the disap- 
pearence of the zero-bias anomalies under anneal- 
ing, by assuming that the presumed static disorder 
anneals away at room temperature. However, this 
suggestion fails a simple quantitative consistency 
check: let us model the constriction region by a Cu 
cylinder 40 nm in diamter and 40 nm long, with 
/ = 3 nm. Estimating the resistance of this cylin- 
der using the Drude model yields R — 7fl, which 
would be the dominant part of the resistance of the 
device {R < lOil in the lower- i? devices). If anneal- 
ing now removes sufficient disorder that the ZBA 
disappears, I would have to increase considerably, 
implying that the overall R of the device would nec- 
essarily decrease by tens of percent, which contra- 
dicts (Cu.3b) (according to which resistance some- 
times even increases under annealing). 

(c) According to (Cu.4b) even a somewhat smaller 
amount of disorder {1^7 nm) than assumed by 
WAM has been observed to cause voltage depen- 
dent conductance fluctuations due to quantum in- 
terference. However, these tell-tale signs of static 



20 



disorder were never seen in the quenched ZBA sam- 
ples (Cu.4c), though they did appear as soon as dis- 
order was purposefully induced using electromigra- 
tion (Cu.4b). In other words, in Cu nanoconstric- 
tions the signature of static disorder is conductance 
fluctuations, not a ZBA. 

3. In the static disorder scenario, the conductance de- 
pends only on the average disorder in the bowl 
(not on the precise configuration of individual de- 
fects). Therefore, it is unclear how to account for 
the complex behavior of the ZBA under thermal 
cycling (Cu.3c), under electromigration (Ti.4), and 
for the F-dependence of the slow ZBA fluctuations 
AG{V) observed in metallic glass break junctions 
(MG.2). Particularly the latter two facts clearly 
demonstrate that the ZBA strongly depends on in- 
dividual defects. 

4. The static disorder scenario predicts a H^^^ be- 
havior for the magnetoconductance, and hence is 
at odds with the very weak H dependence (Ti.5) of 
the Ti ZBAs. 

5. The static disorder scenario provides no hint at all 
about the possible origin of the conductance tran- 
sitions. WAM have suggested that these may be 
due to superconducting regions in the constriction 
(caused by an attractive electron-electron interac- 
tion at short range), but this suggestion fails to 
account for the presence of several different transi- 
tions in the same sample (moreover, superconduci- 
tivity in a Cu sample seems highly implausible). 



2. Magnetic Impurities 

The asymptotic dependence of the conductance on In V 
or InT (Cu.7) is reminiscent of the magnetic Kondo ef- 
fect, where the resistance increases as InT with decreas- 
ing T (as long as T > Tr-). However, there are at least 
three strong arguments that rule out magnetic impurities 
as the source of the anomalies: 

1. An effect due to magnetic impurities would not an- 
neal away at higher temperatures (Cu.3a), since 
magnetic impurities are stable within constrictions, 
not annealing away at room temperatures over a 
time scale of 6 months p2] . 

2. If the magnetic Kondo effect were at work, a mag- 
netic field would cause a well-known Zeeman split- 
ting in the zero-bias conductance dip, as has been 
observed in nanoconstrictions intentionally doped 
with the magnetic impurity Mn [Q section 5.2], as 
shown in Fig. H(a). However, in the devices under 
present consideration, a Zeeman splitting has never 
been observed (Cu.8b). 



Magnetic impurities in metal break junctions have 
been observed to cause ZBAs that do not exhibit 
splitting because the Kondo ternperature scale is 
larger than the Zeeman energy [ p3[ . However the 
ZBA signals caused by these impurities are very 
different than the ones we investigate, because they 
exhibit Fermi liquid scaling (a = 2) rather than the 
a = 0.5 we measure. 



3. TLS Population Spectroscopy 

Many of the qualitative features of the quenched Cu 
ZBAs can be understood within the framework of Kozub 
and Kulik's (KK) theory of TLS-population spectroscopy 
||l^,^, which has recently been extended by Kozub, 
Rudin, and Schober |l^. This theory assumes that the 
constriction contains TLSs (labelled by i) with non-zero 
energy splittings A^ , so that the application of a voltage 
will induce a non-equilibrium population ni±{V) of the 
higher and lower states \±)i of each TLS. Assuming that 
these two states have different cross-sections af for scat- 
tering electrons, the resistance R{V) will then depend 
non-linearly on the voltage and temperture. According 
to Kozub and Kulik, the differential resistance has the 
form 



1 dR 
RdV 



^ 2^-^"^ " ^7) tanh(ir-i)^(t.„ r„ g,) 



(Al) 



where Vj = eV/Aj, tj — k^T/Aj and Cj and qj are geo- 
metrical coefficients depending on the location of the i-th 
TLS in the constriction. The function S{v,T,q), which 
they calculated explicitly, determines the shape of the dif- 
ferential resistance curve (see Fig. 2 of [0), which can 
vary quite significantly, depending on the parameters q 
and T. Note also that since the signs of {u^ — aj) are 
arbitrary (except under special assumptions, see |]64|), 
Eq. ( |A1| ) predicts that ZBAs of both signs should occur. 
The shape of the ZBAs measured by RB is qualita- 
tively the same as that predicted by KK's theory when 
the latter is averaged over many impurities (KK's theory 
has significant freedom for curve-fitting, due to the un- 
determined parameters qj and Aj). However, since the 
function depends on the two parameters v and r sepa- 
rately, this theory cannot account for the existance of a 
scaling law found in (Cu.6a), and certainly not for the 
specific value a = ^ of the scaling exp onent (Cu.6b). 
Moreover, as pointed out in section VH, recent related 
experiments by Keijsers et al. ||2^, in which the behavior 
of individual fast TLSs were probed, shows that they can 
not be reconciled with KK's theory |pfl] (see also pTf). 
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4. Properties of External Circuit 

It was pointed out to us by G. Schon pa that fluctua- 
tions in the voltage V due to fluctuations in the external 
circuit can be shown to lead to a conductance G{V,T) 
that satisfies the V/T-scalmg behavior in (Cu.6) (but 
with the exponent a determined by the external resis- 
tance of the circuit, and hence non- universal) . 

However, since the ZBAs only occur in quenched sam- 
ples (Cu.l), and since ZBAs anneal away at room temper- 
ature (Cu.3a), they must be due to some internal prop- 
erties of the sample. Hence they cannot be due to prop- 
erties of the external electrical circuit, such as external 
voltage fluctations. 



5. Charge Traps and Other Possibilities 

The insulating material used in the devices, namely 
amorphous Si3Ni4, may contain charge traps pq] , which 
could act as Anderson impurities or quantum dots 
through which conduction electrons could hop. This 
could cause dips in the differential conductance through 
several mechanisms, such as Kondo scattering from 
Anderson impurities p7[ |, inelastic hopping conduction 
||,|9| or Coulomb blockade effects [||. 

However, charge traps can be ruled out for the present 
experiments for the following reason. A charge trap has 
in fact been unambiguously observed in a different ex- 
periment by Ralph and Buhrman pq ]. The conduc- 
tance shows a very characteristic peak at 1^ = 0, in 
complete contrast to the ZBA-dip. The suggestion of 
Ref. iQ that this is a Kondo peak that can be asso- 
ciated with Anderson hopping of electrons through the 
trap was taken up by Konig et al. [[71| , who calculated 
the conductance G(V, T) for this scenario and found rea- 
sonably good agreement with that experiment. In other 
words, if charge traps are present, their signals are un- 
mistakable, and very different from the ZBAs of present 
interest. 

Other reasons ruling out charge traps as causes for the 
ZBAs may be found in Ref. |lll. Section 6.6.2. Also in 
Ref. pi[ . Section 6.6, a number of other mechanisms were 
also considered and ruled out as causes for the observed 
ZBAs: electronic surface states or quasi-localized states 
within the metal, defect rearrangement, mechanical in- 
stabilities, superconducting phases and heating effects. 



APPENDIX B: MAGNETIC FIELD 
DEPENDENCE IN 2CK SCENARIO 



which the 2CK scenario could conceivably produce an H- 
dependence for the ZBA. These are the iJ-tuning of the 
asymmetry energy Az {H) , and channel symmetry break- 
ing. Both drive the system away from the degenerate 
2CK fixed point (but not in precisely the same manner), 
so that H enters as a relevant perturbation. However, 
we shall conclude that both mechanism are too weak to 
explain the strength of the observed iJ-dependence in 
quenched Cu samples. 



1. //-Tuning of A 

One possible mechanism by which H could couple to 
the system is by tuning |40 4l|] the TLS asymmetry en- 
ergy Az {H) , and hence the energy splitting A(H) , which 
are then random functions of H (see section V Al ) . With 
A( H) as a relevant perturbation, the analysis of sec- 
tion 



VIC applies directly, and a correction to the con- 



It was stated in section V C 2 that (contrary to the in- 
terpretation we had previously offered in Ref. [||) 2CK 
physics is unable to account for the strong magnetic field 
dependence (Cu.8a) of the ZBA. To illustrate this, we 
now investigate the two most obvious mechanims through 



ductance proportional to A(iJ) can be expected. 

However, in this scenario, the magnetoconductance 
G{H) should be a random function of H (since A(_ff) 
is), whereas it seems to be always positive for the samples 
investigated in more detail. Note also that it would be 
incorrect to attribute the non-universal non-monotonic 
features seen at large H for sample #2 to the random 
behavior of A(if ), since closer scrutiny reveals that this 
behavior is due to the //-motion of the conductance tran- 
sitions (Cu.9b,v). Moreover, since //-tuning of A has its 
origin in quantum interference, it is expected to occur 
mainly in strongly disordered environments, which the 
Cu samples are decidedely not [see (Cu.5)]. Furthermore, 
it would cause conductance changes of order 2e^//i per 
TLS substantially smaller than those observed (Cu.8a) 
(particularly since the signs of the conductance changes 
for different TLSs are random, leading to partial cancel- 
lations). 

Hence, it seems as though //-tuning of A is not con- 
sistent with the observed /J-dependence of the quenched 
Cu samples. 



2. Channel Symmetry Breaking by H 

The second mechanism by which a magnetic field could 
affect a 2CK system is Pauli paramagnetism, which 
breaks channel symmetry (recall that the channel index 
a refers to the Pauli spin |, j.) by causing a net magnetic 
moment M = filHN{£p) [||, eq, (10.11)]. Any such 
symmetry-breaking term can in principle give corrections 
to the critical behavior, and should hence be included in 
the CFT analysis. 

A channel-symmetry breaking field is known to be a 
relevant perturbation with scaling dimension — i (see 
eqs. (3.15) of p^). Hence, in dire ct ana logy to our anal- 
ysis of the effects of A in section VIC , it causes a per- 
turbation (to be denoted by ?7i = 3) to the conductance 
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and 



G{V,T,H), described by Eq. (|l]) with a^ = 

^3 = ^3I^b\H\/ E^ . Here A3 is a dimensionless number 
of order unity, and E3 is an energy that sets the scale at 
which \H\ becomes important. Only the absolute value 
of H enters, because the model is otherwise symmetric 
in spin "f, |, so that the sign of the channel-symmetry- 
breaking field cannot be important. 

This correction term in Eq. (El]) implies that at a fixed, 
small temperature To and V = 0, the conductance obeys 



G{0,To,H)-G{0,To,0)oc\H\ 



(Bl) 



As was argued in Ref. [^, the available data is at least 
qualitatively not in contradiction with this prediction, 
since Fig. |4| shows non-analyticity at H = and an 
initial roughly linear behavior (Cu.8a) (note though, that 
il^/^ behavior can not be ruled out either). 

What are the effects of channel-symmetry breaking at 
sufhciently large if? Presumably, the polarization of the 
Fermi sea will become so strong that one channel of con- 
duction electrons (the one with higher Zeeman energy) 
will decouple froni-the impurity altogether, and the sys- 
tem will cross oveila to the one-channel Kondo fixed point, 
at which the conductance T-exponent is a = 2. Hence, at 
this fixed point the conductance, at fixed, large H, should 
obey the V/T scaling relation Eq. (|I|), with a = 2 |^, 
eq. (D29)]. 

However, it seems unlikely that these considerations of 
the large- iJ regime have any relevance at all for the Cu 
samples, since at large magnetic fields, the conductance 
transitions have moved into the ZBA-regime (Cu.9b,v), 
presumably destroying all remnants of universal 2CK 
physics. (Indeed, a scaling analysis for sample #2 at 
fixed H = 6 T shows best scaling at neither a — ■^ nor 2, 
but at a = 0.3 |g^, though this value probably does not 
have special significance either.) 

Having investigated the phenomenology to be ex- 
pected from channel symmetry breaking, let us now step 
back and estimate the likely magnitude of this effect. 
The Pauli paramagnetism that causes channel symme- 
try breaking merely shifts the bottom of the spin-up 
Fermi sea relative to that of the spin-down Fermi sea 
by fisH, while their respective Fermi-surfaces remain 
aligned. Therefore, the magnitude of the effect that 
H has on the Kondo physics near the Fermi surface 
will be of order figH/D (where D ^ e^ is the band- 
width), and hence negligible. Though this argument is 
not conclusive (e.g. in poor-man scaling approaches the 
band-width is renormalized to much smaller values of or- 
der D' = max[y, T, A]), it casts serious doubts on the 



The cross-over behavior between the fixed points [e.g. the 
behavior of G(0, T, H)] can not be calculated from CFT, 
which can only describe the neighborhood of fixed points; it 
might be possible, though, to calculate this function exactly 
using Bethe-Ansatz techniques. 



channel-symmetry breaking scenario, in particular be- 
cause the observed amplitude of the magnetoconductance 
is by no means small. 

Thus, we have to conclude that 2CK physics cannot 
acco unt fo r the observed ii-behavior (Cu.8a). In sec- 



tion 



VC2| it was therefore suggested to be linked to the 

as 



_ff- motion of the conductance transitions [Vc{H) 
H increases, (Cu.9b,v)]. 
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FIG. 1. Cross-sectional schematic of a metal nanoconstric- 
tion. The hole at the lower edge of the Si3N4 is so small that 
this region completely dominates the resistance of the device. 
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FIG. 2. A typical conductance curve for a constriction con- 
taining structural defects: (a) The upper curve, showing a dip 
in conductance at F = and voltage-symmetric spikes, is the 
differential conductance for an unannealed Cu sample at 4.2 
K. The lower curve shows the conductance of the same device 
at 4.2 K, after annealing at room temperature for 2 days. The 
curves are not artificially offset; annealing changes the over- 
all conductance of the device by less than 0.5%. (b) Point 
contact phonon spectrum at 2 K for the device before anneal 
(dashed line) and after anneal (solid line). 
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FIG. 3. (a) Differential conductance versus voltage at 4.2 K 
for a Cu sample which underwent repeated thermal cycling 
[[ill. The time sequence runs from the bottom curve to the 
top. Curves are artificially offset. The first 2 excursions were 
to 77 K, the next 5 to room temperature, (b) Differential 
conductance for a Cu sample intentionally doped with 6 % 
Au. Static impurities reduce the electronic mean free path but 
completely eliminate the zero-bias anomaly of interest to us. 
(c) Differential conductance and (d) point contact spectrum 
for a Cu device at 1.8 K in which disorder has been created by 
electromigration (which means that a high bias (100-500 mV) 
has been applied at low temperatures so that Cu atoms moved 
around). 



FIG. 4. (a) Conductance signals for 500 ppm magnetic Mn 
impurities in Cu at 100 mK, showing Zeeman splitting in an 
applied magnetic field, (b) The ZBA signals from quenched 
Cu samples exhibit no Zeeman splitting, demonstrating that 
they are not due to a magnetic impurity. However, the shape 
and amplitude of the ZBA does depend on magnetic field. 
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FIG. 6. A generic two-level-system, with (bare) energy 
asymmetry A^ and tunneling rate A^,. An electron- assisted 
tunneling event is depicted: an electron scatters of the TLS 
and induces the atom to tunnel. 
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FIG. 5. (a) V-dependence of the differential conductance 
for B — and T — 100 mK, for three different samples #1, 
#2, 7^3. (b) T-dependence of the conductance for 5 = and 
V — for the three samples of (a), and a fourth. 
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FIG. 7. (a) Temperature dependence of the V — conduc- 
tance [G(0, T) — G(0, To)] for the 4 unannealed Cu samples 
of Fig. H plotted versus T^/^ The values of G(0, T) for the 
different samples, extrapolated to T = as shown are for sam- 
ple #1: 2829 e^/h, sample #2: 3973 e^/h, sample #3: 30.8 
e^//i, and sample #4: approximately 2810 e^ /h. (b) Voltage 
dependence of the differential conductance at T = 100 mK for 
some of the same samples as in (a), plotted versus V^'^ . The 
size of deviations from T^ " behavior in (a) (1 part in 3000) is 
consistent with the magnitude of amplifier drift in these mea- 
surements, as they were performed over several days. The 
T/-dependent measurements in (b) are less subjective to such 
drift problems, as they are taken over a much shorter time 
span. 



FIG. 8. (a) Voltage dependence of the differential conduc- 
tance for sample #1 of Fig. M, plotted for temperatures rang- 
ing from 100 mK (bottom curve) to 5.7 K (top curve), (b) 
The same data, rescaled according to Eq. (Ilj) with a = i, 
and plotted vs. v^'^ — (eV/fcar)^". The low-temperature, 
low-voltage data collapse onto a single curve [linear for large 



,1/2 



in agreement with Eq. (E3|)], with deviations when the 



voltage exceeds 1 mV. 



28 



0.2 



0.2 



0.5 
a 




0.5 
a 



I , 




-5.0 



-2.5 0.0 2.5 5.0 

sqrt(eV/kBT) 




0.0 5.0 

sqrt(eV/kBT) 



FIG. 9. The deviation parameter D{a) of Eq. (pi), which 
quantifies the quality of scaling, for (a) sample #1 and (b) 
sample #2. The minimum of D{a) defines the value of a that 
gives the best scaling, giving a = 0.48 ± 0.05 for sample #1 
and a = 0.52 ± 0.05 for sample #2. 



FIG. 10. (a) Differential conductance data for sample #2, 
at temperatures from 200 mK to 5.7 K, and (b) for sam- 
ple #3, at temperatures from 50 mK to 7.6 K, rescaled ac- 
cording to Eq. and plotted vs. v^^'^ = (eV/fcijT)^''^. 
The low-voltage, low-temperature data collapse well onto one 
curve for sample ^2, but not for sample #3, partly due to the 
existence of TLSs with Kondo temperatures within (rather 
than above) the temperature range of the measurement. 
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FIG. 11. Representative conductance curves which lie 
along the scaling curves for samples #1, #2 and #3 of each 
of the three samples in Figs. H(b), Fig. IWa) and (b), respec- 
tively. For sample #1, the curve corresponds to T=l.l K, 
for sample ^2 1.4 K, and for sample #3 250 mK. The rea- 
son for selecting these particular curves was that among all 
those lying along the scaling curve, they had the best sig- 
nal-to-noise ratio for each sample, (a) The j/-axis is scaled by 
the value of B^ determined from the temperature dependence 
of the V = conductance for each sample (values listed in Ta- 
ble 6.1). (b) In addition, the a;-axis is scaled with a number 
fli for each sample. The scaling curves for all three samples 
seem to lie on one universal curve. 



FIG. 12. (a) Transition voltage Vc{H) of the conductance 
transition |32]as a function of magnetic field for a quenched 
Cu constriction at 4.2 K, showing bifurcation (Cu.9b,iv). At 
high fields Vc{H) — > (Cu.9b,v), the dependence on H being 
quadratic, (b) Vc{H) for five other samples at 4.2 K, with 
associated decreases to 0, quadratically in H. 
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FIG. 13. The DC conductance Goc (as opposed to dif- 
ferential conductance G used in all other plots) at several 
magnetic fields for a quenched Cu sample [B2| . An increasing 
magnetic field broadens the conductance transition at Vc and 
moves Vc toward zero voltage, destroying the enhancement at 
V = Vc of Gdc above the normal conductance. Note that 
although a large applied magnetic field eliminates the con- 
ductance transition, a V = minimum in the conductance 
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FIG. 14. Magnetic field dependence of the V=0 conduc- 
tance for the 3 unannealed Cu samples at 100 mK. (a) Ab- 
solute magnetoconductance. (b) Magnetoconductance scaled 
by the value of Bs for each sample, (c) Magnetoconductance 
relative to the change in conductance between 100 mK and 
6 K. An applied magnetic field alters, but does not eliminate, 
the zero-bias conductance signal due to TLSs. 
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FIG. 15. The differential conductance for a 20 f2 Ti 
constriction (sample 4), for temperatures of 6.0, 4.0, 2.0 
and 1.4K, scales well when plotted in the scaling form 
[G{V,T) - G{0,T)]/T^^^ vs. {eV/kBT)'^^^ of Eq. ^) with 
o ~ ^ [see (Ti.2)]. Top inset: the unsealed conductance 
AG(y) = G{V,T)-G{0, 1.4K). Bottom inset: The deviation 
parameter D{a) of Eq. (|24[), which quantifies the quality of 
scaling, according to which a — 0.52 ± 0.05. 
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FIG. 16. (a) AG{V) = G{V,T) - G(0,0.76mK) vs. 
V^^'^ for a 19 n Ti constriction (sample 5) at 1.41K and 
76mK. (b) The same data do not scale when plotted as 
[G{V,T) -G(0,T)]/Ti/2 vs. {eV/kBTy^^, due to "prema- 
ture saturation" (see arrow) of the 76mK curve, as discussed 
in (Ti.3). 
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FIG. 18. The magnetic field dependence (at and 5T) at 
lOOmK for two Ti constrictfons [see (Ti.5)]: (a) sample 5; (b) 
sample 6. Data for (a) are offset by e^ /h at V = for clarity. 



FIG. 17. (a) The ZBA of the Ti sample 5 before and after 
electromigration [see (Ti.4)], which changed the saturation 
energy T^ of Eq. (H) from 2.3 to 1.4K. (b) A 22 f7 Ti con- 
striction (sample 6) shows two distinct saturation energies 
Tj,! and Tx2, which change upon electromigration from 0.9 to 
1.5K and 9.7 to 6.8K, respectively. 
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TABLE I. Measured parameters of the scaling functions for 
the ZBAs in 3 Cu constrictions (samples 1 to 3) and one Ti 



constriction (sample 4). B-e 
and Fi is dimensionless. 



Fo and _F\ have units K e /h, 



# 


5e 


Fo 


Fi 


^1 = sir 


1 


7.8 ±0.2 


4.2 ±0.3 


-5.7 ±0.9 


-0.73 ±0.11 


2 


25.2 ±0.7 


12.8 ±0.8 


-19.7 ±1.5 


-0.78 ±0.06 


3 


10.3 ±0.4 


6.0 ±0.6 


-7.7 ±1.6 


-0.75 ±0.16 


4 


1.69 ±0.2 


0.89 ±0.05 


-1.37 ±0.1 


-0.81 ±0.10 



0.75 




8 10 



FIG. 19. ZBA fluctuations AG(V) [see (MG.2)] due to the 
presence of a slow fluctuator in two metallic glass constric- 
tions studied in Ref. po|. The squares give AG{V) for Fig. 2, 
curve 3 of po|; the triangles give the noise amplitude mul- 
tiplied by 2 (for visibility) of Fig. 4, curve 1 of [EOHuncer- 
tainties ~ 0.1e^//i). The flts were obtained in Ref. p(\\ using 
VZ's theory to calculate AG{V). For the solid curve a single 
fast TLS was assumed, experiencing modulations in asym- 
metry energy between A^ = 8 meV and A^ = 3 meV, with 
Tk = 17 K. For the dashed curve two fast TLSs were as- 
sumed, with A^i = 9,Aii = 6.2 meV, Tm = 8.9 K, and 
Arf = 4.2, A',2 = 2.8 meV, Tk2 = 6.2 K. 
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